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Abstract. Universal bi-Hamiltonian hierarchies of group-invariant (multi- 
component) soliton equations are derived from non-stretching geometric curve 
flows j(t, x) in Riemannian symmetric spaces M = G/H, including compact 
semisimple Lie groups M = K for G = K X K, H = diag G. The derivation of 
these soliton hierarchies utilizes a moving parallel frame and connection 1-form 
along the curve flows, related to the Klein geometry of the Lie group G D H 
where H is the local frame structure group. The soliton equations arise in 
explicit form from the induced flow on the frame components of the principal 
normal vector N = V x ~y x along each curve, and display invariance under the 
equivalence subgroup in H that preserves the unit tangent vector T = ■jx in the 
framing at any point a; on a curve. Their bi-Hamiltonian intcgrability struc- 
ture is shown to be geometrically encoded in the Cartan structure equations 
for torsion and curvature of the parallel frame and its connection 1-form in the 
tangent space T 7 M of the curve flow. The hierarchies include group-invariant 
versions of sine-Gordon (SG) and modified Korteweg-de Vries (mKdV) soliton 
equations that are found to be universally given by curve flows describing non- 
stretching wave maps and mKdV analogs of non-strctching Schrodinger maps 
on G/H. These results provide a geometric interpretation and explicit bi- 
Hamiltonian formulation for many known multicomponent soliton equations. 
Moreover, all examples of group-invariant (multicomponent) soliton equations 
given by the present geometric framework can be constructed in an explicit 
fashion based on Cartan's classification of symmetric spaces. 



1. Introduction and Overview 

The theory of integrable soliton equations displays many deep links to differ- 
ential geometry, particularly as found in the study of geometric curve flows. In 
this paper, group-invariant soliton equations and their bi-Hamiltonian integrabil- 
ity structure are derived from studying non-stretching flows of curves in symmet- 
ric spaces G/H. Such spaces describe curved G-invariant Riemannian manifolds, 
generalizing the classical two-dimensional Riemannian geometries i.e., the sphere 
S 2 ~ SO{3)/SO(2), the hyperbolic plane H 2 ~ SL(2, R)/SO(2), and the Euclidean 
plane K 2 ~ Euc(2)/SO(2), which are characterized by constant (positive, negative, 
zero) curvature. 

Geometric curve flows in constant curvature spaces M = S 2 ,H 2 ,M. 2 are well 
known to yield [14, 16, 22, 18] the mKdV hierarchy of scalar soliton equations and 
its hereditary recursion operator. The starting point is to formulate the flow of a 
non-stretching (inextensible) curve j(t, x) as geometrically given by jt — h±N+h\\T 
in an adapted moving frame (T, N) = (7^, *~/ x ) (i.e. a Frenet frame) along the curve, 
where x is arclength and * is the Hodge-star operator on vectors in the tangent 
plane M 2 ~ T X M. Preservation of the non-stretching condition 7^ • r y x = 1 under the 
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curve flow yields the relation d x h\\ = K,h±, allowing the tangential component h\\ of 
the flow to be expressed in terms of the normal component h± and the curvature 
invariant n — N ■ V^T of the curve. The flow equation combined with the Serret- 
Frenet structure equations of the moving frame then determine an evolution of this 
curvature invariant Kt = lZ(h±) + \h± where 1Z = D + n 2 + k x D x 1 k is seen to be 
the hereditary recursion operator of the mKdV hierarchy of soliton equations [29] , 
and x is the Gaussian curvature (+1,-1,0) of M. In particular, the flow given by 
h±_ = k x , = i« 2 yields the mKdV evolution equation on k to within a convective 
term, n t — x K x = K3x + \k 2 k x (called the +1 flow in the hierarchy). Higher order 
mKdV evolution equations arise from linear combinations of flows h± = TZ u (k x ), 

= D^ 1 (nlZ n (k x )) , n = 1, 2, . . .. The entire mKdV hierarchy n t — x K x — 
thus sits in the class of geometric flows in which the evolution of the curve in the 
normal direction is a function hj_ = /ij_(/t, n x , . . .) of the differential invariants of the 
curve [25, 26]. Such flows correspond to geometric equations -f t = f(j x , V X 7 X , . . .) 
satisfied by the curve 7(^,2;). 

It is less widely known that the geometric curve equations produced in this 
manner from the mKdV hierarchy of evolutions on k are close analogs of Schrodingcr 
maps 7t = J(7x), J — iV^, on S 2 ~ CP 1 (identifying i with * after complexification 
of the tangent plane). For instance the mKdV evolution itself yields -f t — k x N + 
\k 2 T = V x (k *~f x ) + ^k 2 j x describing [1] a geometric curve equation j t — ^7(7^)7 
K 1 = V x + ||V a; 7 a ;| 2 as obtained through the Serret-Frenet equations V^T = kN , 
V^iV = — kT. This flow operator K 1 is the mKdV analog of a modified Schrodingcr 
operator J 7 = iV^, + {Qxg^ x ) x that preserves |7 X | = 1. One geometrical difference 
between these operators is that 7^ • J 7 (7 X ) = whereas 7^ • Kj^x) = i« 2 > 0, 
implying (by dimensional considerations) that the non-stretching Schrodingcr flow 
of the curve -){t,x) in S 2 is actually stationary -f t — J 7 (7 2 ) = (so thus K t = 0), 
while in contrast the non-stretching mKdV flow j t — K 7 ("/ x ) is dynamical. 

A related geometric evolution [1] comes from the kernel of the mKdV recursion 
operator lZ(h±) = giving the evolution equation n t = x^± (called the —1 flow 
in the mKdV hierarchy), with h± satisfying d x h± = — nh\\, d x h\\ = nh±. This 
flow obeys the conservation law = d x (h 2 + h\) = V x (7t • 7t) which implies 
it is conformally equivalent (under rescalings of t) to a flow with uniform speed 
|7i| = \Jh 2 + h\ = 1. The evolution of the curve is then given by h± = sin(8), h\\ = 

cos(9) in terms of the nonlocal invariant 9 = — J ndx of 7, yielding the sine-Gordon 
equation 9 tx — — xsin(0) induced from the evolution on k [21, 28]. Equivalently, 
the relations 8 = arcsin(K t /x) and 9 X = —k yield a hyperbolic equation Kt x — 
— K\/x 2 ~ K t 0n K - The corresponding geometric curve equation satisfied by j(t, x) 
is V :E 7t = V x (sin(9)N) + V x (cos(9)T) = after the Serret-Frenet equations are 
used. This curve flow is recognized [1] to be a uniform speed, non-stretching wave 
map V x 7t = V t 7j; = on S 2 , with t, x viewed as light cone coordinates for the 
wave operator. 

Moreover, the geometry of the two-dimensional surfaces swept out by all these 
non-stretching curve flows j(t, x) in M = S 2 ,H 2 ,M. 2 turns out to encode the bi- 
Hamiltonian integrability structure [30] of the mKdV soliton hierarchy. The rel- 
evant geometry [1] is given by the Ricmannian connections V^V^, in the flow 
direction j t and tangent direction j x , specifically that they have vanishing torsion 
^ t lx - V 2 ,7t = [7t,7x] = and carry constant curvature [V t , VJ = R(j t ,Jx) = 
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x(lt • *7a;)* where R(X, Y) is the Riemann curvature tensor of M with X, Y in 
the tangent plane. A projection of these equations into the normal space of 7 
yields n t = Ti(w) + xhj_, zu = expressed in terms of the flow invariant 

zu = N ■ V t T, where H — D x , J = D x + k,D~ x n are found to be respectively 
the mKdV Hamiltonian cosymplectic and symplectic operators for the hierarchy of 
mKdV evolutions on k. In this setting h±d/dK represents a Hamiltonian vector 
field, and zndn represents a variational covector field with w = SH/Sk holding for 
some Hamiltonian H — H(k, k x , . . .). Compatibility of the operators Tt, J implies 
that 1Z — TLoJ and its adjoint TV = JoTi generate a hierarchy of commuting vector 
fields = TZ u {k x ) and corresponding involutive covector fields = TV n {n), 

n = 0,1,2,..., related by hj n) = H(m (n) ), m {n+1) = J(h ± (n) ). For n > 1 
this hierarchy possesses a bi-Hamiltonian structure h±^ = H(5H^ /Sk) = 
J^ 1 {5H ljl+1 ^ 1 1 '5k) where 3~ x is the formal inverse of J defined on the rc-jet 
space of k, and with the local Hamiltonians = H^ n \n, k x , . . .) determined 

by td^ = SH^/Sk. Thus Ti^J^ 1 are formally a bi-Hamiltonian pair of cosym- 
plectic operators. Alternatively, £ — HJTL — D x + k 2 D x + k x D x 1 (kD x ) provides 
an explicit cosymplectic operator compatible with H, giving the bi-Hamiltonian 
structure /i ± (n) ' = £{5H {n - 1] / 5k) = H(5H (n) /5k). 

The bottom of the mKdV hierarchy hj_^ = k x , zu^ = k, = ^k 2 orig- 

inates geometrically from the x-translation vector field K x d/dn which is a sym- 
metry of the operators TL,J. In terms of the jet space variables (x, k, n x , . . .), 
the entire hierarchy has the form of homogeneous polynomials with non-zero scal- 
ing weight under the mKdV scaling symmetry x — > Ax, k — > A _1 k. As a 
consequence, through a scaling formula (cf. [2]) the Hamiltonians are given by 
H {n) = j^(k + xn x )w {n) = jl^D^inh^) = j^h^ (modulo total x- 

derivatives) for n > 0, where d x h^ — Kh± n \ Thus J hiS^dx represents a 
(scaled) Hamiltonian functional on the re-jet space of k. The same expression car- 
ries over to the —1 flow in the mKdV hierarchy, h±^~^ = sin(#), h\\- = cos(0), 
given in terms of k = —9 X . Specifically, hj_^~^ = H(5H^~^ /5k) holds for the 
Hamiltonian H^ 1] = with zu^ 1] = 5H { ~ 1] /5k = D' 1 (5H^ 1} /SO). 

Recent work [3] has generalized these results to the two known vector ver- 
sions of the mKdV equation [34, 6], by deriving their bi-Hamiltonian integrablity 
structure along with their associated hierarchies of (higher order) vector soli- 
ton equations from geometric curve flows in the Riemannian symmetric spaces 
SO(N + l)/SO(N) ~ S n and SU(N)/SO(N). These two spaces exhaust (cf. [19]) 
all examples of irreducible symmetric spaces G/SO(N) given by compact simple Lie 
groups G, and in the case N = 2 they both coincide with the classical Riemannian 
spherical geometry S 2 ~ SO(3)/SO(2) ~ SU(2)/SO{2). 

This geometric derivation of the vector mKdV hierarchies involves several main 
ideas. Firstly, for a non-stretching curve 7, the components of the principal normal 
in an adapted moving frame along 7 play the role of a natural Hamiltonian flow 
variable. Secondly, the torsion and curvature equations of the Riemannian connec- 
tions V t , V on the two-dimensional surface of any flow j(t,x) of such curves in 
G/SO(N) geometrically encode Hamiltonian cosymplectic and symplectic opera- 
tors Tt, J, where x is the arclength on 7. In particular, the Hamiltonian structure 



4 



STEPHEN C. ANCO 



looks simplest if the moving frame along 7 is chosen so that it has a parallel connec- 
tion matrix [9] given by an algebraic version of the basic geometric property V^iV 
|| T, V X T 1 T of an adapted moving frame generalized from S 2 to G/SO(N). 
For such a moving parallel frame the torsion and curvature equations reduce to 
a simple form which can be derived directly from the infinitesimal Klein geome- 
try of the tangent spaces q/so(N) in G/SO(N), where g is the Lie algebra of G. 
Thirdly within the class of non-stretching curves whose parallel moving frame is 
preserved by an 0(N — 1) isotropy subgroup of the local frame structure group 
SO(N) of the frame bundle of G/SO{N), the encoded operators H,£ = HJH 
comprise an 0(N — l)-invariant bi-Hamiltonian pair. Moreover, they possess the 
mKdV scaling symmetry and exhibit symmetry invariance under x-translations. 
As a result, these operators yield a hereditary recursion operator 1Z = H.J gener- 
ating a hierarchy of integrable curve flows in G/SO(N), organized by their scaling 
weight, in which the flow variable satisfies an 0(N — l)-invariant vector evolution 
equation. The +1 flow in the respective hierarchies for G = SO(N + 1),SU(N) 
is given by the two known vector mKdV soliton equations [34], and there is a —1 
flow described by vector hyperbolic equations [6] (variants of vector SG equations) 
coming from the kernel of the respective mKdV recursion operators in the two hier- 
archies. All the flows in both hierarchies correspond to commuting bi-Hamiltonian 
vector fields in the x-jet space of the flow variable. Finally, the evolution equations 
of the curve "f(t,x) produced by the ±1 flows in each hierarchy are G-invariant 
generalizations of the geometric map equations found for these flows in the case 
S 2 ~ SO{3)/SO(2) ~ SU(2)/SO(2), identified as wave maps and mKdV analogs 
of Schrodingcr maps, on G/SO(N). 

More recently there has been an extension of these results [4] deriving complex 
generalizations of the vector mKdV hierarchies along with their bi-Hamiltonian in- 
tegrability structure from geometric curve flows in the Lie groups G = SO(N + 
1),SU(N). The derivation adapts the moving parallel frame formulation of non- 
stretching curve flows to these Lie groups by viewing them as Riemannian symmet- 
ric spaces in the standard manner [20, 19], with the frame invariance group being 
a natural isotropy subgroup U(N — 1) C G contained in the local structure group 
of the frame bundle of G = SO(N + 1), SU(N). This leads to bi-Hamiltonian 
hierarchies of curve flows j(t, x) in which the flow variable (again described by 
the components of the principal normal vector along the curves) obeys a U (N — 1)- 
invariant vector evolution equation. In particular, the hierarchies contain the known 
complex generalizations [34] of the two versions of vector mKdV soliton equations, 
as well as complex generalizations of the corresponding vector SG equations which 
were first obtained by symmetry-integrability classifications [6]. 

A full generalization of these ideas and results will be presented here to obtain bi- 
Hamiltonian hierarchies of group- invariant soliton equations arising from geometric 
curve flows in general Riemannian symmetric spaces G/H. All irreducible examples 
of these spaces divide into two types, where either G is a simple Lie group (and H 
is a compact subgroup invariant under an involutive automorphism of G) or G is a 
Lie group product K x K (and H is a diagonal subgroup) such that G/H ~ K is 
a compact simple Lie group. 

Theorem 1.1: For each Riemannian symmetric space G/H there is a family of 
bi-Hamiltonian hierarchies of non- stretching curve flows described by G-invariant 
geometric maps on G/H. The components of the principal normal vector along 
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these curves in a moving parallel frame satisfy group-invariant (multicomponent) 
soliton equations of mKdV and SG type, with an explicit bi-Hamiltonian structure. 
In each hierarchy the geometric maps corresponding to the SG flow and the lowest- 
order mKdV flow are universally given by a non-stretching wave map and a mKdV 
analog of a non- stretching Schrodinger map. 

Remark 1.2: The invariance group of the soliton equations and their bi- 
Hamiltonian integrability structure consists of the linear isotropy subgroup in 
Ad(H) that preserves the tangent vector for the corresponding curves in the tan- 
gent space T Q G/H = g/\) at the origin o. (By G-invariance, any geometric curve 
flow is equivalent to one that passes through the origin o in G/H.) 

Remark 1.3: Up to isomorphism, the distinct bi-Hamiltonian hierarchies in 
the family admitted by a given Ricmannian symmetric space G/H are in one-to- 
one correspondence with equivalence classes of unit-norm elements contained in any 
fixed maximally abelian subspace in jj/f) modulo the Weyl group W C Ad(H). 

These results provide a geometric origin and unifying interpretation that en- 
compasses many examples of multicomponent soliton equations. All examples can, 
moreover, be written down in explicit form, using Cartan's classification of sym- 
metric spaces [19]. 

To begin, some preliminaries are stated in section 2 on the relevant differential 
geometric and Lie algebraic properties of symmetric spaces [19, 20], particularly 
the infinitesimal Klein geometry [33] attached to g/f), and the Cartan subspaces 
in jj/f) determined by the Lie algebras t),g of H, G. These properties are used in 
an essential way in section 3 for the study of non-stretching curve flows 7(i, x) in 
G/H, where x is arclength (so \^ x \ = 1 is the non-stretching property using the 
Riemannian metric on G/H). The main starting point is given by the pull back of 
the torsion and curvature equations of the Ricmannian connection on G/H to the 
two-dimensional surface of the flow of 7, as derived directly by a frame formulation 
in terms of the infinitesimal Klein geometry of g/f). Based on the decomposition 
of g relative to ad(e) where e is a unit vector belonging to a Cartan subspace in 
g/i), a natural algebraic construction is given for a parallel moving frame along 7, 
evolving under the flow. Then in section 4, the moving frame components of the 
torsion and curvature are shown to encode a bi-Hamiltonian pair of cosymplectic 
and symplectic operators H, J with respect to a Hamiltonian flow variable given 
by the moving frame components of the principal normal along 7. The proof of this 
main result will involve extending the standard theory of bi-Hamiltonian structures 
[29, 15] to the setting of Lie-algebra valued variables. 

The hierarchy of bi-Hamiltonian commuting flows produced by the operators 
TL,J\s studied in section 5, where the +1 and —1 flows are constructed and found to 
yield group-invariant multicomponent versions of mKdV and SG equations associ- 
ated to all irreducible Riemannian symmetric spaces. The geometric map equations 
produced by these universal ±1 flows are derived in detail. Concluding remarks on 
extensions of these results will be made in section 6. 

Some related results of interest have been obtained in recent literature. In one 
direction, Hamiltonian operators for one of the vector mKdV equations as well as for 
a related scalar-vector KdV system have been derived from non-stretching curve 
flows in iV-dimensional constant curvature Ricmannian geometries (i.e. the TV- 
sphere) and in ./V-dimensional flat conformal geometries (i.e. the Mobius iV-sphere), 
using two different approaches. One approach [30, 31] viewed these geometries as 
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locally modeled by Klein geometries on R defined respectively by a Euclidean 
isomctry group action Euc(N) and a Mobius (conformal) isometry group action 
Mob(N) ~ SO(N + 1,1). A subsequent derivation [1, 5] used just the intrinsic 
Riemannian and conformal connections defined on the underlying manifold S N . 
Both derivations utilize the approach of studying curve flows via parallel moving 
frames, first introduced in [22, 27, 30]. The same approach has also been applied 
recently [36] to curve flows in symplectic geometry Sp(N+l)/Sp(N) x Sp(l), giving 
a geometric derivation of a quaternionic (non-commutative) mKdV equation. 

Earlier fundamental work appeared in [7, 8] on multicomponent mKdV equations 
derived from symmetric spaces by algebraic and geometric considerations, and on 
investigation of Hamiltonian structures of such equations. 

Another direction [10, 11] has concentrated on deriving the known scalar soli- 
ton equations from non-stretching curve flows in various classical plane geometries: 
hyperbolic plane, affine and fully-affine plane, plane similarity geometry. These 
geometries together with the Euclidean plane are characterized by their respec- 
tive isometry groups SL(2), SA(2) and A(2), Sim(2), Euc(2), each acting locally 
and effectively on R 2 ; as such they comprise the main examples of planar Klein 
geometries G/H ~ R 2 where G is the isometry group and H C G is a stabilizer 
subgroup that leaves fixed the origin in R 2 . In this geometric setting there is a 
natural group-invariant notion of geometric curve flows based on the evolution of 
differential invariants [29] of a non-stretching curve 7, formulated using a moving 
frame (T,N) = (■y x ,lxx) where a; is a natural G-invariant arclength on 7. General- 
izations to non-stretching flows of space curves in Klein geometries on R 3 have also 
been studied in a few cases [12, 13]: Euclidean space Euc(3), similarity geometry 
Sim(3), affine space A(3), centro-affine space SL(3). The results obtained to-date 
have provided an elegant geometric origin for many scalar soliton equations. For 
instance, the KdV equation and the Sawada-Kotera equation are noted to arise re- 
spectively as a centro-affine and affine version of the mKdV equation coming from 
Euclidean geometry. It should be possible to derive the bi-Hamiltonian integra- 
bility structure of these solitons equations by adapting the parallel moving frame 
approach used in the present work to the setting of affine geometry and similarity 
geometry. 

2. Preliminaries 

A symmetric space M = G/H is irreducible if it is not a product of smaller sym- 
metric spaces. All irreducible Riemannian symmetric spaces have a classification 
into two types [20, 19]: 

(I) G is a connected simple Lie group and H C G is a compact Lie subgroup 
invariant under an involutive automorphism a of G (in this case M is compact or 
noncompact according to whether G is). 

(II) G is a product K x K for a connected compact simple Lie group K and H is 
its diagonal product subgroup in G on which a permutation of the product factors 
acts as an involutive automorphism a, with G/H ~ K ~ H (in this case, since G 
is compact, so is M). 

Henceforth let M = G/H be an irreducible Riemannian symmetric space in 
which G acts effectively on the manifold M by left multiplication. The automor- 
phism a induces a decomposition on the Lie algebra of G, given by a direct sum of 
vector spaces g = [} © m such that cr(t)) = f), <r(m) = — m, with Lie bracket relations 
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[f), t)] g C t), [f),tn] fl Cm, [tn, m] g C f), where I) is identified with the Lie algebra of H 
and m = T M is identified with the tangent space at the origin o in M. The group 
action of G on M provides a canonical isomorphism T X M ~ G m. 

2.1. Riemannian metric and connection. The Riemannian structure of G/H 
comes from the Cartan-Killing form (-, •) of g. Since g is a simple Lie algebra, 
the subspaces f),m are orthogonal (f),m) = and hence the Cartan-Killing form 
restricts to yield a nondegenerate inner product on m which is Ad g (H) invariant. 
This inner product (•, -) m has a definite sign, negative or positive, according to 
whether g is compact or noncompact. In the noncompact case, the Lie algebra 
g = \) ® m is dual to g* = f) ©im which defines a compact real Lie algebra contained 
in the complexification of g. This leads to a corresponding duality [20, 19] between 
compact and noncompact Riemannian symmetric spaces of type I. As there are no 
noncompact Riemannian symmetric spaces of type II, we will henceforth restrict 
attention to the compact case in order to include both types I and II. 

Under the canonical identification m ~ G T X M, the negative-definite inner prod- 
uct (•, -) m defines a Riemannian metric tensor g on M given by 

g(X,Y) = -(X,Y) m (1) 

for all G- invariant vector fields X, Y in T X M . Associated to this metric g is a unique 
torsion- free G- invariant Riemannian connection V, determined by Vg = 0, where 

T(X,Y) = V X Y-V Y X- [X,Y] = (2) 

is the torsion-free property. Here [•, •] denotes the commutator of vector fields on 
M. The Riemann curvature tensor of the connection V is given by 

R(X,Y) = [V X ,V Y ] - V [X<Y] - -ad m ([A,y] m ) (3) 

with ad m viewed as a linear map on T X M ~ G m. This curvature tensor is G- 
invariant and covariantly constant Vi? = 0. Note that left multiplication by G 
acting on M generates the isometries of the manifold M, namely the isomctry 
group of g is G. 

Recall, the linear isotropy group H* of G/H is the representation of H by linear 
transformations on T M at the origin o in M . This group is isomorphic to Ad m (H) 
through the identification T M = m. 

2.2. Klein geometry and soldering frames. For any Riemannian symmetric 
space M = G/H, the Lie group G has the structure of a principal iJ-bundle over M, 
tt : G — > G/H [20, 33]. Consider the Maurer-Cartan form wg on G, which is a left- 
invariant g- valued 1-form that provides the canonical identification T C G = g where e 
is the identity element of G. The Maurer-Cartan form satisfies dwe + \ [wg, wg] b — 
and hence, viewed geometrically, it is a zero-curvature connection on the bundle 
G. This collective structure is called [33] the Klein geometry of M = G/H. It 
has a useful reformulation locally on M involving the Lie algebra structure of the 
tangent spaces of G/H. 

Let ip : U — > 7r _1 (C7) C G be any local section of this bundle, where {U,ip) 
are coordinates for a local trivialization G w U x H. The tangent space of the 
trivialization U x H at any point (x,ip(x)) is isomorphic to the Lie algebra g of G. 
In particular, note tt* :g~mxl)^mas vector spaces, with T X M ~ G m = g/f), 
T,p( x )H ~ f) and T( x ,ip(x))G — = m [). Now, the pullback of the Maurer-Cartan 
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form wq by ip yields a g-valued 1-form b lo (locally) on M, satisfying the Cartan 
structure equation 

^+|[VH = o (4) 

where D is the total exterior derivative operator on M and [■, -] B is understood to 
denote the Lie algebra bracket in q composed with the wedge product in T*M. 
(Informally, we can view B io : T X M ~ G m = g/t) — ► g as attaching an infinitesimal 
Klein geometry to g/f).) Note a change of the local section tp, described by a smooth 
function h : U — > H yielding the local section jp = tph, induces a change of b lo given 

by 

b lo = Adih-^cj + h^Dh. (5) 

This is called [33] a gauge transformation on b uj relative to the gauge tp. 
Introduce 

e:=i(« W -a(M). u> ~\(>u> + v(>u>)), (6) 

denoting the even and odd parts of s u>, defined using the automorphism a of q. 
The corresponding decomposition under a gauge transformation on B ui yields 

e = Ad(h~ 1 )e, Q = Adt/r 1 )^ + hr x Dh. (7) 

Proposition 2.1: OnU C M, the m-valued 1-form e represents a linear coframe 
and the t)-valued 1-form lo represents a linear connection, with H as the structure 
(gauge) group defining a bundle of linear frames. 

Let the 2-forms 

% = De + [u,e] g = Ve and $H = Du + = \D,V] (8) 

denote the m-valued torsion and rj-valued curvature associated to the linear con- 
nection and coframe, where 

V = D+[w,\ (9) 

is the (gauge) covariant exterior derivative. 

Proposition 2.2: The Cartan structure equation (4) decomposes with respect 
to g to give 

1 = 0, $H = -i[ e ,e] m , (10) 

so thus T> is torsion-free and has covariantly- constant curvature. 

The precise interpretation of e and lo in these propositions comes from a basis 
expansion e = e a m a , lj = u^h 1 , where m a ,h l are any fixed orthonormal basis for 
m, f) with respect to the Cartan-Killing inner product, 

(m Q ,m h ) = -5 ab , (h*,h J ') = -5 ij , (m a ,h 4 )=0. (11) 

Then {e a } will define frame covectors spanning the cotangent space of M at x, and 
{LJj} will define connection 1-forms related to the exterior derivative of the frame 
covectors by the torsion-free property of the covariant derivative V on M at x. 

Most importantly, e and lo provide a natural soldering of this frame geometry 
of M = G/H locally onto the Riemannian geometry of M = G/H. Let e* be dual 
to e, namely a m-valued vector such that — (e*,e) m = id^ is the identity map on 
T X M. Then e* represents an orthonormal linear frame on U C M, so thus {e*} 
as given by a basis expansion e* = e*m a will define frame vectors spanning the 
tangent space of M at x and obeying the orthonormality property e*Je& = 6 a b- 
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Theorem 2.3: The linear coframe e and linear connection u> determine a Rie- 
mannian metric and Riemannian connection on M = G/H given by 

g = -(e,e) m , V- = -(e*,V(e\-)) m = -<e*,D(eJ-) + [w, (ej-)] fl )m. (12) 

In particular, for all X, Y in T X M , 

g(X,Y) = -{e x ,e Y ) m (13) 

and 

V X Y = — (e*,V x e Y ) m , or equivalently T> x e Y = e\\7 x Y (14) 
where e x = e\X,e y = e\Y are the coframe projections of X,Y. Moreover, the 
Riemannian curvature is determined by the curvature of the exterior covariant de- 
rivative V, 

R(X, Y)Z = -(e*, [V x ,V Y ]e z - V [XY] e z ) m = -(e*, [K XY , e z ] s ) m (15) 
with (cf Proposition 2.2) 

m xY =m\(X AY) = -[e x ,e Y ] m . (16) 

Hence 

R(X,Y) = ([[e x ,e Y ] m ,e] g ,e*) m , (17) 

or equivalently 

#(•,•) = -<[ej-,ej-] m ,[e*,e] m >„. (18) 

Here note [e*,e] m defines a f)-valued linear map on T X M displaying skew sym- 
metry with respect to g: g(Y, [e*, e] m X) = [e Y , e x ] = -g(X, [e*,e] m Y). 

Remark: The torsion-free property of the Riemannian connection is expressed 
by (cf Proposition 2.2) 

T(X,Y) = -(e*, {V Y e)\X- (V x e)\Y) m = -(e",Z x , Y ) m = (19) 

with 

% X , Y =%\(XAY). (20) 
The soldering structure of the Riemannian connection has the following alterna- 
tive formulation. 

Corollary 2.4: For all X, Y in T X M , 

V x e=[e,w x ] fl , [V x ,V Y ]e-V [XY] e=-e\R(X,Y) = [e,m x ^] s (21) 

where uo x = X\uj is the linear connection projected in the X direction on M. 

In examples, it will be convenient to work in a matrix representation for m, 
viewed as subspaces in gl(N, R) for suitable N > 1. Hence m a , h 1 will then represent 
N x N matrices satisfying the commutator relations 

[h\ h 3 } = c t: > k h k , [m a , m b ] = c ah M l , [h l , m a ] = c la b m b , (22) 

where c^k, c ah iC ta b are the structure constants of q with respect to the matrix basis 
{m a ,h 2 }. More directly, () and g can be represented as matrix subalgebras in 
gl{N, R), with m represented as the matrix subspacc g/f) in gl(N, M). Then the 
soldering structure is given by matrix formulas 

g(X,Y) = - X tv(e x e Y ), (23) 
V X Y = - X tr(e*V x e Y ), R(X, Y) = - X tv([e* ,e][e x ,e Y }), (24) 

for some constant \ depending on the choice of matrix representation of m, where 
e x = (e a \X)m a , e Y = (e a jy")m a arc N x N matrices belonging to m; e = e a m a , 
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e* = e*m a are N x N vector-valued matrices belonging to T*M ® m, T X M ® m 
respectively. Moreover, 

V x < - -c^^^et, (25) 

R(X, Y)e* a = -c a a (m x Y ) ie l = c cd lC lb a X c Y d el (26) 

with X a = e a \X,Y a = e a \Y. 

Due to the soldering, the frame covectors obey the orthonormality property 
gK^tr 1 = -tr(m a m b ). 

2.3. Action of ad(e). It is useful to consider an algebraic decomposition of g with 
respect to ad(e) where e is the soldering frame. Let X be a vector in T X M and 
consider the element e x = e\X given by the frame projection of X into m C g. 
The centralizcr of e x is the set c(e x ) of all elements annihilated by ad(e x ) in the 
Lie algebra g, namely [e x ,c(e x )] = 0. This set has the obvious structure of a 
Lie subalgebra of g, since ad(e x )[c(e x ), c(e x )] = from the fact that ad acts as 
a derivation on the Lie bracket of g. Write c(e x ) m , c(e x )[, for the intersection of 
c(e x ) with the subspaces m, \) in g = m© fj, and c(e x )^, c(e x )j|- for their respective 
orthogonal complements (perp spaces) where 

m = c(e x ) m ®c(e x )i, f) = c(e x )„ 8 c(e x )^. (27) 

Thus, note [e x ,c(e x ) m ] = [e x ,c(e x )(,] = 0, while [e x ,c(e x )^] ^ 0, [e x ,c(e x )(,] 7^ 
0. 

Proposition 2.5: 77ie centralizer subspaces c(e x ) m , c(e x )(, o/e x ftoue tfte Lie 
bracket relations 



[c(e x )[,,c(e x ) b ] C c(e x )(,, 
[c(e x )i,,c(e x ) m ] C c(e x ) m , 
Hence the perp spaces satisfy the relations 
[c(e x )„,c(e x )^] C c(e x )^-, 
[c(e x )„,c(e x )4] C c(e x )^, 



(28) 

[c(e x ) m , c(e x ) m ] C c(e x V (29) 

[c(e x ) m ,c(e x )^]Cc(e x )^ (30) 
[c(e x ) m , c(e x )j}-] C c(e x )^. (31) 



The proof uses just the Lie bracket structure of g along with the invariance 
of the Cartan-Killing inner product. For instance, consider [c(e x )(,, c(e x )f,]. Be- 
cause ad(e x ) is a derivation it annihilates [c(e x )(,, c(e x )(,] and hence, since f) is a 
Lie subalgebra of g, [c(e x )(,, c(e x )(,] is contained in c(e x )(,. Thus [c(e x )f,, c(e x )f,] 
is orthogonal to c(e x )jj-, and so by invariance, = (c(e x )jj-, [c(e x )(,, c(e x )f,])f, = 
(c(e x )[,,[c(e x )(,,c(e x )^-]>[, which implies [c(e x )f,, c(e x )^] is contained in c(e x )j}-. 
The other relations are established similarly. 

The observation that e x lies in c(e x ) m , combined with the Lie bracket relations 
in Proposition 2.5, shows that ad(e x ) maps c(e x )f, into c(e x ) m , c(e x )^ into c(e x )4, 
and vice versa. Consequently, ad(e x ) 2 is well-dchncd as a linear mapping of the 
subspaces c(e x )jj- and c(e x )^ into themselves. These subspaces are likewise taken 
into themselves by all the linear maps ad(c(e x )f,), corresponding to the infinitesimal 
action of the group of linear transformations that preserves e x . In particular, let 
H* x be the closed subgroup in the linear isotropy group H* given by Ad(ft,)e x = e x . 
Since [ad(e x ), ad(c(e x )f,)] = ad([e x , c(e x )(,] = 0, note that c(e x )(, is isomorphic to 
the Lie algebra of H* x and that the linear transformations in this group commute 
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with the linear map ad(e x ) 2 . Schur's lemma applied to ad(e x ) 2 now gives the 
following result. 

Proposition 2.6: (i) c(e x )^ and c(e x )^ are isomorphic as vector spaces under 
the linear map a,d(e x ). (ii) c(e x )^ and c(e x )^ each decompose into a direct sum of 
subspaces given by irreducible representations of the group H* x on which the linear 
map ad{e x ) 2 is a multiple of the identity. 

Note dimc(e x )^ = dimc(e x )jj- = dimm — dm\H* x , so these spaces will be of 
minimal dimension (depending on e x ) precisely when the group H* x is of maximal 
size. 

The decomposition of these vector spaces in Proposition 2.6 will be maxi- 
mally irreducible whenever the linear map ad(e x ) 2 has equal eigenvalues, whereby 
ad(e x ) 2 — xidj^ is a multiple of the identity map id± on c(e x )^ ~ c(e x )^. In 
this situation, additional Lie bracket relations can be established as follows. In- 
troduce the bracket (•,•) := [ad(e x )- , • ] on the centralizer subspace c(e x )g := 

c ( e x)m ® c i e x)i) °f e x m S' an d write Ix ■= ad(e x ) 2 . Since ad(e x ) acts as a 
derivation, note 

Ix(y,z) = [a,d(e x )Ixy,z] + [a,d(e x )y,I x z]+2[I x y,a,d(e x )z] 

= {Ixy,z) + (y,I x z)-2{z,I x y) (32) 

for all vectors y, z in c(e x )g. Now suppose Ix = xid±, where x is a non-zero 
constant. Then the identity (32) yields xid^(y, z) — 2\{y,z) — 2x(z,y) and hence 
(y,z)± = 2(z,y)± where a subscript _L denotes projection onto c(e x )g. This iden- 
tity implies {y,y)± = 0, and hence by polarization, (y,z)± = 0. As a result, it 
follows that = (y, z)± - (z, y)± = [ad(e x )y, z]± + [y, ad(e x )z}± = id ± ad(e x )[y, z] 
whence [y, z]± — since ad(e x ) is nondegenerate on c(e x )^. 

Proposition 2.7: If ad(e x ) 2 = \id^ is multiple of the identity on c(e x )„ ~ 
c ( e x)t) ' Ihcn these centralizer subspaces obey the Lie bracket relations 

[<e x )^, c(e x )£] C c(e x )„, [e(e x )i, c{e x )i\ C c(e x )„, (33) 
[c(e x )i,c(e x )^Cc(e x ) m . (34) 



2.4. Cartan spaces. The algebraic structure of the spaces c(e x )^ and c(e x )^ 
along with the group H* x will have a key role in the later construction of a parallel 
moving frame formulation for curve flows in Riemannian symmetric spaces. Part 
of this algebraic structure is effectively independent of e x . 

Let co(e x ) denote the center of c(e x ), which is an abelian subalgebra of Q. 
Write Co{e x ) m , Co{e x )f,, for the intersections of Co(e x ) with m, [), and write a for 
any maximal abelian subspace in m. 

It can be shown [19] that the image of any maximal abelian subspace a under 
the group H* is the vector space Ad(H)a = m. Hence e x is contained in some such 
subspace, namely a = c (e x ) m C c(e x ) m C m. These subspaces a C m are known 
to have the following characterization [19]. 

Lemma 2.8: (i) Any two maximal abelian subspaces a are isomorphic to one 
another under some linear transformation in H* ~ Ad(H). (ii) Every maximal 
abelian subspace a is invariant under the Weyl group W — C(a)/N(a), which is 
a finite subgroup contained in the linear isotropy group H* , where C(a) and N(a) 
are respectively the centralizer group and normalizer group of a in H* . (Hi) Any 
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maximal abelian subspace a is the — 1 eigenspace in some Gartan subalgebra of the 
Lie algebra g = m © \] under the grading cr(m) = — m, er(f)) = f), and hence a is the 
centralizer of its elements mm, a = c(a) C m. 

Henceforth the subspaces a C m will be referred to as Gartan subspaces of m in 
g. The close relationship between the subspaces o and the Cartan subalgebras of g 
leads to a more explicit description of the algebraic decomposition of g with respect 
to &d(e x ) 2 . 

Let g c be the complexification of g. Recall, a root of g is a linear function, 
a, from any Cartan subalgebra g° of g c into C, such that for all vectors w in g°, 
a,d(w) has an eigenvector with the eigenvalue ce(w). Let the set of all roots of g be 
denoted A. The following lemma is adapted from results in [19], utilizing the fact 
that c(e x )^ is contained in a- 1 as a consequence of the inclusion a C c(e x ) m . 

Lemma 2.9: Let A m be the set of roots a of g such that a(a) ^ for some a in 
a Cartan subspace a C m. Since g is compact, a(a) will be an imaginary function 
(namely a : a — > iM.). Then there exists a basis of vectors {y a } in c(e x )^ and 
{z a } in t(e x )ff, indexed by a certain subset of the roots a in A m , satisfying the Lie 
bracket relations 

[ya,a] g = \a{a)z a , [z a ,a] g = -ia(a)y a (35) 

for all vectors a in a. Hence the eigenvalues of the linear map ad(a) 2 on its 
eigenspaces in c(e x ) 1 | l , c(e x )j|- are given by a(a) 2 < (cf Proposition 2.6). Note 
ad(a) 2 will be a multiple of the identity precisely when its eigenvalues a(a) 2 are 
equal. 



3. Parallel moving frames and bi-Hamiltonian operators 

On a n-dimensional Ricmannian manifold M, a frame is a basis of the tangent 
space T X M at each point x in M, and a coframe is a dual basis of T*M . Let 
{e a } be an orthonormal coframe and {e*} its dual frame, a — l,...,n, namely 
g{e a ,eb) = e*Jef, = 5 a b- Covariant derivatives of the frame in any direction X at 
a point in M are given in terms of the frame connection 1-forms oj a b = —uJb a by 
V x e a = — (u>a b \X)ei,. Here the frame structure group (modulo reflections) is the 
rotation group SO(n). A moving frame along a curve 7(2) in M is the restriction of 
a frame to T 7 M, or correspondingly a coframe to T*M, describing a cross section 
of the orthonormal frame bundle over 7. SO(n) gauge transformations allow a 
moving frame to be adapted to 7 such that 

{(e±)a, a = 2,...,n 

where = X is the unit tangent vector along 7 and each (e*j_) a is a unit normal 
vector with respect to 7. An adapted frame is preserved by the SO(n — 1) group 
of local rotations in the normal space of 7. This gauge freedom can be used to 
geometrically adapt the connection matrix 

(cj x ) ab :=u ab ]X = g(e* a ,V x e* b ) (37) 

to 7 such that the derivative of the tangent vector in the frame is normal to 7 while 
the derivative of each normal vector is tangential to 7. Such a moving frame is said 
to be parallel]^)]: 

J_ e|, V x (el) a || ef. (38) 
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Note the connection matrix u> x of a parallel moving frame is skew-symmetric with 
only its top row and left column being nonzero, as a consequence of the parallel 
property 

2(4V x ep = 0, 5((el) a ,V x (el) fe ) = 0. (39) 

This property is preserved by any rigid (x- independent) SO (n — 1) rotation in the 
normal space of 7, which thereby defines the equivalence group for parallel moving 
frames. 

Any parallel moving frame differs from a classical Frenct frame [17], if n > 2, 
by an x-dependent SO(n — 1) rotation acting in the normal space of the curve 
7 (see [35] for a geometric interpretation of this rotation). It is well known that 
the components of a Frenet connection matrix along 7 arc differential invariants of 
the curve [17]. Consequently, the components of a parallel connection matrix are 
thus differential covariants of 7 in sense that they are invariantly defined up to the 
covariant action of the equivalence group of rigid SO(n — 1) rotations. 

There is a purely algebraic characterization of the connection matrix of a parallel 
moving frame [3]. Fix the standard orthonormal basis of R™, given by n row 
vectors m a whose o th component is 1, with all other components equal to 0, where 
a = 1, . . . , n. Similarly fix the standard basis of the Lie algebra so(n) induced from 
the isomorphism R™ Al n ~ so(n), where a wedge denotes the vector outer product. 
(Note the natural action of so(n) on R™ is defined here by right multiplication of 
skew matrices on row vectors.) This basis consists of n(n — l)/2 skew matrices 
h* = m° A m b whose only nonzero components are 1 in the a th row and 6 th column, 
and —1 in the 6 th row and a th column, where i is identified with (a, b) such that 
a < b, for a,b — 1, . . . , n. 

Now, there are associated to the frame e* and its covariant derivatives the R n - 
valued linear coframe e = e a m a and the so(n)-valued linear connection 1-form 
uj = uj a bm a A m b satisfying soldering relations [20] analogous to the ones (12), 
where e represents a linear map from R™ onto T X M. These relations state that V 
has no torsion 

de-eAw = I=0 (40) 

and carries curvature 

dw+[w,w]„ (n) =9t (41) 
determined by the frame components of the Ricmann curvature tensor [V,V] of 
the metric g on M. 

Write rri|| :— m 1 when a = 1, m" := m a when a = 2, . . . , n, so thus e\\ = eim^ 
and e± = e a m" describe the tangential and normal parts of the linear coframe 
e = e || + ej_ relative to X, with ejjX = for the tangent vector X and eiijY = 
for all normal vectors Y along 7. Consider the derivatives of the coframe e along 7 
and decompose the connection matrix uj x := lo\X = u> x + ui x into separate parts 
given by uj° x = (w^J^m" A = and u x = (u! x )i a m^ A m" 7^ 0. Observe 
that in this decomposition := m" A is a basis for the Lie subalgebra 
so(n - 1) ~ R" _1 A M™ _1 where M™ -1 c M™ is the subspace orthogonal to rri||, 
while hj| := rri|| A m" is a basis for the perp space of so(n — 1) C so(n). 

Proposition 3.1: A linear coframe e = e^ + e± along a curve 7(x) in a Rie- 
mannian manifold M is a parallel moving frame iff its connection matrix lu x van- 
ishes on the Lie subalgebra so(n — 1) C So(n) of the SO(n—l) rotation subgroup that 
leaves the tangent covector e\\ invariant under the SO{n) frame structure group. 
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This algebraic characterization of the connection matrix of a moving parallel 
frame for curves in Ricmannian geometry has a direct generalization to the setting 
of Klein geometry following the ideas in [3] . 

3.1. Construction of parallel moving frames in Klein geometry. Let 7(x) 

be a smooth curve in some n-dimensional compact Riemannian symmetric space 
M = G/H, and write X — -y x for its tangent vector with an afhne parametrization 
so that x is the arclength. Locally on M , introduce a m-valued linear coframe e 
and a F)-valued linear connection u> satisfying the soldering relations (12) in terms 
of the Riemannian metric g and connection V on the manifold M, where m = g/t) 
is the vector space quotient of the Lie algebras g, t). This structure naturally arises 
from the Klein geometry of the Lie group G viewed as a principle bundle over M 
whose structure group is the Lie subgroup H C G. The projections of e and u> 
along the curve 7(2;) 

e x :=e\X, w x :=w\X (42) 

provide a moving frame formulation of this curve such that the frame structure 
group is the linear isotropy group H* ~ Ad(if) acting as local gauge transfor- 
mations (7). The arclength property of x in this formulation is expressed by 
\X\ g = \e x \ m = 1 so e x has unit norm in the Cartan-Killing metric on m (where 
l e xlm = — ( e x> e jsf)m)- Informally e can be regarded as a nondegenerate cross 
section of the Klein bundle (T*M ® m, H) over 7(2;), where m ~ R n . This bun- 
dle is a reduction [20] of the orthonormal frame bundle of M, with the structure 
group SO(n) reduced to H. (In fact, H C SO(n) with equality holding only when 
M = G/H is a constant curvature manifold, M = S n ~ SO(n + l)/SO(n) in the 
compact case.) 

Recall [19], the vector space m = g/f) has the algebraic structure of a union of 
orbits generated by the group H* acting on any fixed Cartan space (a maximal 
abelian subspace) a c m (all of which are isomorphic). Every such subspace a, 
moreover, is a finite union of closed convex sets a*(m) C a that are mutually 
isomorphic under the Weyl group W C H* leaving a invariant. (Up to closure, 
these sets 0* in o each can be identified with Weyl chambers [19] associated to the 
set of positive roots of the Lie algebra g relative to m.) 

Since H* acts as a group of local gauge transformations on e and u, a gauge 
equivalent coframe exists such that, at every point x along 7, e x is a constant 
unit vector that lies in the (fixed) Cartan space a C m on which the Weyl group 
acts as a discrete gauge equivalence group. This leads to the following algebraic 
characterization for curves 7(2;). 

Proposition 3.2: The constant unit vector e x C a defines an algebraic invari- 
ant of the curve 7. Thus, all curves on M = G/H can be divided into algebraic 
classes in one-to-one correspondence with the unit-norm elements in the closed con- 
vex set a*(m) ~ a/W as determined by e x . 

For a given vector e x , a linear coframe e along 7 is determined only up to local 
gauge transformations (7) given by the group H* x C H* that preserves e x C a. 
Such transformations will be referred to as the isotropy gauge group of 7. The 
isotropy gauge freedom can be used to adapt e and w to 7 in a preferred way, gen- 
eralizing the algebraic notion of a parallel moving frame for curves in Riemannian 
geometry stated in Proposition 3.1. 

Definition 3.3: A m-valued linear coframe e along a curve 7(2;) in a Riemann- 
ian symmetric space M = G/H is said to be H-parallel if its associated [)-valued 
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linear connection lu x in the tangent direction X — j x of the curve lies in the 
perp space c(e x )^ of the Lie subalgebra c(e x )^ C i) of the isotropy gauge group 
H* x C H* that preserves e x given by a fixed unit vector in o C m. 

To explain the geometrical meaning of such a coframe, first decompose the dual 
frame e* = ej| + and the connection w = + lu 1 - according to 

m = c(e x ) m ®c{e x )i, I) = c(e^ c(e^, (43) 

relative to the unit vector e x C a C m. The m-valued linear frame e* will provide a 
framing of the curve 7 in M (namely a linear map from the vector space m* ~ K™ 
onto T*M over 7(0;)) once a basis expansion for m, f) is introduced. Note the 
condition for e* to be ii-parallel is 

lu x := J!jX = (44) 

along with 

D x e x =0, \e x \ m = l. (45) 

Remark: This property u> x = is maintained by rigid (x-independent) 
gauge transformations in H* x that preserve e x , since under these transformations, 
Ad(h~ 1 )X = X implies Ad(/i _1 )cjj f belongs to c{e x \ and hence the transformed 
linear connection (7) satisfies uj x = 0. All such transformations comprise the 
equivalence group, denoted Ht, for i?-parallel moving frames. The infinitesimal 
equivalence group is simply given by the centralizer Lie subalgebra c(e x )^ C t), 
as seen from the Lie bracket relations (31) and (30). Note fl| preserves the de- 
composition (43) of the tangent space T 1 M ~ G — m into algebraically parallel and 
perpendicular parts with respect to X . 

Let {mjj, hy} be respective orthonormal bases for the subspaces c(e x ) m , c(e x )^, 

so 

ad(e x )mjj = ad(e x )hj| = (46) 

(where a = 1, . . . , dimc(e x ) m ; i — 1, . . . , dim c(e x )%). Recall, the perp subspaces 
c ( e v)mi c i e x)\) are isomorphic under the linear map &d(e x ), and so let {m" ,h^} 

be respective orthonormal bases for c(e x )^ ~ c(e x )^ related by = ad(e x ) _1 hY 
(where a,i = 1, . . . , dim c(e x )^ = dim c(e x )jj-). Orthonormality of each basis is 
given by 

(mf,m^) m = (ml,mi) m = -<5 afc , (47) 

(h^hj^HhiX)^-^, (48) 

and 

(mi,m||) m = (hi,hj) f ,=0. (49) 

Now, write 

(e|)„ = -<mjf, e*) m , (el)„ = -<mi, e*) m , (50) 

and 

{J\) i = -{bi ] ,u) f) , (w- L ) < = -<hi,w) b) (51) 

yielding the basis expansion of e*,u. As a consequence of properties (46)-(49), the 
vectors (ej|) a and (e* ± ) a will be said, respectively, to be algebraically parallel (||) and 
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perpendicular (_L) to the tangent vector X = ^ x along the curve. This framing is 
adapted to 7 if it contains X as one of the vectors parallel to X, (say for a = 1) 

(e|)i=X, (52) 

which can be achieved by letting mjj = e x be one of the basis vectors of c(e x ) m . 

Lemma 3.4: The frame vectors {(e|) , (e^) a } provide an orthonormal basis 
of TjM over the curve j(x). Through the soldering formula (25), their covari- 
ant derivatives in the tangent direction X — j x are determined by the connection 
matrices (uJ x )iC th a with (u x )iC lb a = 0. 

The Lie bracket relation (31) between c(e x )jj- and c(e x ) m now leads to an explicit 
geometric interpretation of the TJ-parallel property. Note Y is algebraically || to X 
iff &d(e x )e Y — [e x ,e Y ] — namely e Y := e\Y belongs to c(e x ) m . Likewise Y is 
algebraically _L to X iff e Y belongs to c(e x )^. 

Theorem 3.5: For any H -parallel frame along a curve j(x), the covariant 
derivatives V x (ep a , a = 1, . . . , dim c(e x ) m , are algebraically _L to the tangent vector 
X = -f x . The other covariant derivatives V x (e^) a , a = 1, . . . , dim c(e x )^, are 
algebraically || to X iff the Lie bracket relation [c(e x )jjj, c(e x )jj-] C c(e x ) m holds, 
which is always the case whenever the linear map ad(e^) 2 is a multiple of the 
identity on the vector spaces c(e x )^ ~ c(e x )^- (cf Proposition 2.4)- 

These geometric properties of an f/-parallel frame are a strict generalization 
of the Riemannian case (38) whenever dimc(e x )^ > 1, or [c(e x )^,c(e x )^] g n 
<e x )i ± 0. 

There is a further geometric meaning for an _ff-parallel connection itself. Recall, 
the principal normal along any afBne-parametrized curve is the normal vector [17] 

- V xlx . (53) 

Proposition 3.6: In an H-parallel adapted frame along a curve j(x), the 
principal normal is given by N = (&d(e x )uJ x , e*) m with the frame components 
N a = e a \N = (oJ x )iC ll a := (u> x ) a where c 4l a = ad(e x ) a 4 is an invertible ma- 
trix (a,i = 1, . . . , dim c(e x )^ = dimc(e x )^J. Moreover, these components are 
differential covariants of the curve with respect to the rigid equivalence group ( of 
x-independent linear transformations) HT, C H* (preserving the framing). Namely 

the set of components (w^)iC jl a is invariantly defined by 7(2) up to the covariant 
action of the group H*^ . 

An important final remark now is that the if-parallel condition (44) on the 
linear connection u> can always be achieved by means of a suitable gauge trans- 
formation (7) as follows. Under the isotropy gauge group H* x , uj x is transformed 
to uj x = Ad(/i _1 )w|| f + h~ 1 D x h where h is an arbitrary smooth function from the 
curve 7(2;) C M into the Lie subgroup of H* x C H, satisfying Ad(h~ 1 )e x = 0. 
Condition (44) applied to u x yields an ODE on h(x), h~ 1 D x h = -Ad(/i _1 )w|| f , 
for which local existence of a solution is a standard result. (In particular, in the 
adjoint representation of H and f), this ODE becomes a linear homogeneous matrix 
differential equation D x Ad(h) = — Ad(h)ad(u> x ).) This result is a straightforward 
generalization of the standard existence of a parallel moving frame in the Riemann- 
ian setting. 
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Proposition 3.7: By means of a gauge transformation in the isotropy gauge 
group H* x , an iJ-parallel frame {(e|) , (e^) a } exists for any curve j(x) in a given 
Riemannian symmetric space G/H and is unique up to the equivalence group C 
H* consisting of rigid (x-independent) gauge transformations that preserve the 
algebraic structure (43)-(45) of the framing along 7(2;). 

4. Derivation of bi-Hamiltonian operators from non-stretching 

curve flows 

Consider an arbitrary smooth flow j(t, x) of any smooth curve in an n- 
dimcnsional Riemannian symmetric space M = G/H. Write X — ^ x for the 
tangent vector and Y — -f t for the evolution vector along the curve. Provided the 
flow is transverse to the curve, namely X and Y are not parallel vectors, then -f(t, x) 
will describe a smooth two-dimensional surface immersed in the manifold M. 

The torsion and curvature equations of the Riemannian connection V on M 
restricted to this surface *y(t, x) are given by 

T{X,Y) = V xlt -V tlx = Q, (54) 
R{X,Y) = [V x , V t ] = -ad ra ([ 7 x,7t]m), (55) 

where [•, -] m is the Lie bracket on m = g/f) with the canonical identification m ~ G 
T X M. There is a straightforward frame formulation for these equations obtained 
from the pullback of the Cartan structure equations (10) for any m-valued linear 
coframe e and [)-valued linear connection u) on the surface ^(t, x): 

D x e t ~ D t e x + k,et] 8 - [u t , e x ] g = 0, (56) 

D x u> t - D t uj x + [w x ,w t ]f, = -[e x ,e t ] m , (57) 

where 

e x := e\X = e\^ x , e t := e\Y = ej7 t , (58) 

oj x := uj\X = uj\~f x , uj t := oj\Y = u\^ t . (59) 

These structure equations coincide with the zero-curvature equation (4) for G 
viewed as a i7-bundle in a local trivialization over the surface -f(t, x) in M. In 
particular, in local coordinates for G w U x H, the g-valucd connection 1-form 
e + uj = s uj coincides with the pullback of the left invariant Maurer-Cartan form on 
G to C/| 7 C M, locally giving a soldering of the frame geometry of G/H onto the 
Riemannian geometry of the manifold M. 

4.1. Non-stretching geometric curve flows. A curve flow j{t,x) is non- 
stretching (inextensible) if it preserves the arclength ds = \-f x \dx at every point 
on the curve, namely Z^^l = or equivalcntly V t "/ X (= V ' x 7t) is orthogonal to 
"f x in the Riemannian metric g on M. Now let 7(i, x) be a non-stretching curve 
flow given in terms of an affine parametrization \ j x \ = 1, without loss of generality, 
where x measures arclength along the curve. For any such flow, there exists an 
ii-parallel frame formulation obtained through Proposition 3.7. In particular, a 
suitable local gauge transformation in H* ~ Ad(H) (given by a change of coordi- 
nates for the local trivialization of G) can be applied to the linear coframe e and 
linear connection lo to produce an ii-parallel moving frame along each curve in the 
flow 7(t, x), depending smoothly on t. 
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Lemma 4.1: In an H -parallel frame formulation for non- stretching curve flows, 
the projections (58) and (59) of e and u in the tangent direction X = j x and the 
flow direction Y — j t have the following properties: 

(i) e x belongs to some Cartan space a C m 7 is constant (D x e x — 0), and has the 
normalization \e x \ m = 1; 

(ii) e t := hu + h± and u) t := ro" + w 1 - each decompose into algebraically parallel 
and perpendicular parts relative to e x , namely fty and ro" belong to the centralizer 
subspaces rri|| :— c(e x ) m and hp := c(e x )j, while h± and w 1 - belong to the perp 
spaces mj_ := c(e x )^ and f)j_ := c{e x )^; 

(Hi) u :— lo x belongs to the perp space f)j_; 

(iv) h±, e x , u respectively determine the perpendicular flow vector Y± — (7t)i = 
— (h±,e*) m , the tangent vector X = j x = — (e x ,e*) m , and the principal normal 
vector N = V ' x -f x = — (ad(e x )u, e*) m for the curve. 

The Lie bracket relations in Proposition 2.5 on the spaces m = mil mi and 
t) = f) ii ©f)j_ lead to a corresponding decomposition of the Cartan structure equations 
(56) and (57) expressed in terms of the variables e x , h\\ , h±,w^ , ru ± ,u. The torsion 
equation becomes 

= D x h l{ +[u,h A _] l] , (60) 

= ad(e :c )w ± + D x h± + [u, h\\] + [u, h±]±, (61) 

while the curvature equation becomes 

= D x m'k + [u,m ± ] ||, (62) 

ad(e x )h± = D t u- D x m^ + [zu^u] + [txj^,u]±, (63) 

where [-, -Jy denotes the restriction of the Lie bracket on g to m|| or fjy, and likewise 
[•, denotes the restriction to m_i_ or hj_. Note that, in the present notation, the 
Lie bracket relations consist of 

[m||,m||]Cf,| h ft||,m||]Cm| h [f),,, F)„] C (64) 

[tn||,mj_] C f)j_, [f)||,mj_] C m_L, [l)±,m\\] C mj_, [{jy, f)±] C (65) 

so thus 

[m||,m||]_L = [f) M , m.||]_L = fty, fj M ] _i_ = 0, (66) 

and 

[m,|,m_L]|| = [{)||,mi.]|| = [fu,m||]|| = [f)||,fu]|| = 0. (67) 

Therefore the only brackets with a possibly nontrivial decomposition are [mi , mj_] 
and [f)j_, t)i[. 

Lemma 4.2: T/ie Cartan structure equations for any H -parallel linear coframe e 
and linear connection u pulled back to the two-dimensional surface -f(t, x) determine 
a flow on u = u) x given by 

u t = D x w^ + [u, tu"] + [u, w ± ]± + &d(e x )h±, (68) 
w 1 - = -adi^)- 1 (D x h± + [u,h\\] + [u,h±]±), (69) 

with 

h^-D-^h^, a 7ll = -DJ 1 [«, n 7- L ]||, (70) 
where D" 1 denotes the formal inverse of the total x-derivative operator D . 
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Note that, once h±_ is specified, this determines both the flow on u and the curve 
flow of 7 itself, where h\\ , vu^ , zu 1 - are given in terms of h± by equations (69)-(70). 
A natural class of geometric flows will now be introduced. 

Definition 4.3: A curve flow is said to be geometric if h± is an equivariant 
m^-valued function of (x,u,u x , . . .) under the action of the rigid (cc-independent) 
group of linear transformations HZ , where the arclength variable x is an invariant 
of the curve, and the components of flow variable u are differential covariants of 
the curve (cf Proposition 3.6). 

In the case of polynomial geometric flows, m^-equi variance means that the func- 
tion h± is constructiblc in terms of (u,u x , . . .) together with e x and x using just 
the Lie bracket and Cartan-Killing inner product on these variables. Such flows 
correspond to curve flows j(t,x) such that 74 = f(x,~/ x ,~/ xx ,...) is constructed 
using the metric g(-, •) and the linear map ad x (-) 2 on T X M ~g m via the relations 
(i) to (iv) in Lemma 4.1. 

Proposition 4.4: For a given polynomial geometric flow on u, the curve 7 obeys 
a G-invariant flow equation in M determined through the identifications 

lx <-» e x , V x <-> D x + [u, -} m := V x , ad x (•)<-» ad m (-), g{-, ■) <-> -(•, -) m . (71) 



4.2. Hamiltonian operators and bi-Hamiltonian structure. As a main re- 
sult, the flow equation on u will now be shown to possess an elegant bi-Hamiltonian 
structure. To begin, recall that &d(e x ) as a linear map gives an isomorphism be- 
tween the spaces and fy±. Write h 1 - := &d(e x )h± for the image of h± which 
belongs to the same perp space f)j_ as u. Expressing u in terms of h- 1 and zu 1 - then 
leads to the following operator form for the flow equation (68): 

u t =H(vu ± ) + h ± , vu ± =J(h ± ), (72) 

where 

H = lC\ h± , J = -ad( e J- 1 /C| m ,ad(eJ- 1 (73) 

are linear operators which act on f)j_ -valued functions and are invariant under Htt, 
as defined in terms of the operator 

IC:=D x + [u r } ± -[u,D- 1 [u,-} ll }. (74) 

Theorem 4.5: Ti, J are compatible Hamiltonian cosymplectic and symplectic 
operators with respect to the t) ^-valued flow variable u. For geometric flows given by 
h ± (x, u, u x , . . .), the flow equation on u has an H*^-invariant form u t = Wi( 1 7(/i ± )) 
where Hi = H + J^ 1 is a formal Hamiltonian operator. Consequently, such a 
geometric flow will have a Hamiltonian structure 

u t =Hi{Sf)/5u) (75) 

for some Hamiltonian functional Sj = J H(x, u,u x , . . .)dx iff 

SSj/Su = zu ± =J{h- L ). (76) 

Moreover, all such Hamiltonian flows have a second Hamiltonian structure 



u t = H 2 (S3/6u) 



(77) 
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given by the ( compatible) formal Hamiltonian operator H.2 = TZ T~ti, with a Hamil- 
tonian functional 3 = f F(x, u, u x , . . .)dx such that 

S3/Su = K*(m ± ) = J(K(h^)), (78) 

where 1Z = T-LJ is a (hereditary) recursion operator (andlZ* = JTL is its adjoint). 

Here SS) /Su denotes the variational derivative of the real- valued functional Sj. 
More precisely, 59) /Su is defined as the unique f)j_ -valued function of (x, u, u x , . . .) 
such that 

r (Sf)/Su,h J -)dx = S h ±Sj (79) 



formally holds for all f)j_-valued functions h,- 1 , where S h ± — pr(/i J - • d/du) denotes 
the variation induced by h^- , namely the prolongation of the given vector field 
h 1 - ■ d/du on J°° , and where (-,■) denotes the (negative-definite) Cartan-Killing 
inner product on t)±. (A dot will denote summation over Lie-algebra components 
with respect to any fixed basis.) 

The definition of cosymplectic and symplectic operators and a proof of this 
theorem will be given in sections 4.3-4.5. A summary of the basic theory of bi- 
Hamiltonian operators and Hamiltonian structures for scalar evolution equations 
is given in [29, 15]. The presentation here will extend this theory to a Lie-algebra 
valued setting. 

One main property characterizing Hamiltonian operators consists of the sym- 
plectic structure they induce on the x-jet space (x,u,u x , . . .) := J°° of the flow 
variable u. In particular, there is associated to H the Poisson bracket defined by 

{Sj, <5} n = - J (H(S®/Su), SS)/Su)dx (80) 

for any pair of functionals fj, 6. This bracket is skew-symmetric 

= (81) 

and obeys the Jacobi identity 

{Si, {0, %} h }h + {©, {3, K}h}h + {5, {Si, ®}n}n = (82) 

as a consequence of H being cosymplectic. Let h ± (x, u, u x , . . .) be a function with 
values in the space ty±. Then, on the x-jet space J°°, h ■ d/du is a Hamiltonian 
vector field with respect to Ti if there exists a (Hamiltonian) functional S) such that 

s h ±e = {e,si} n (83) 

for all functionals (3. Associated to J is a symplectic 2-form defined by 



u(hi ■ d/du, ■ d/du) j = - J {hi, 



J{h^))dx (84) 



for any pair of Hamiltonian vector fields given by hi ,h% ■ The property that u> is 
symplectic corresponds to the skew-symmetry and Jacobi identity of the Poisson 
bracket {-,-}j-i defined using the formal inverse of J. (Consequently, J^ 1 is 
formally cosymplectic, and conversely Hr x is formally symplectic.) This 2-form 
<jj gives a canonical pairing between Hamiltonian vector fields and corresponding 
covector fields defined as follows. Let tn- L (x, u, u x , . . .) be a function with values in 
the space i)±. Then the covector field vj 1 - ■ du is dual to a Hamiltonian vector field 
h 1 - ■ d/du if 

w ± -du:=u(- ,h^-d/du)j (85) 
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holds with respect to J. As a consequence of J being symplectic, such covector 
fields are variational, namely there exists a (Hamiltonian) functional £ such that 

5 h ±<£ = -J{m^,h^)dx (86) 

for all vector fields h 1 - ■ d/du. 

Proposition 4.6: The Hamiltonian (co symplectic and symplectic) operators 
(73) give mappings h 1 - ■ d/du — ► J(h ) ■ du and w ■ du — » Hfa- 1 ) ■ d/du between 
Hamiltonian vector fields and variational covector fields on the x-jet space of the 
flow variable u. 

Theorem 4.5 and Proposition 4.6 provide a broad generalization of known results 
on the geometric origin of bi-Hamiltonian structures in Ricmannian symmetric 
spaces M = G/H. Special cases were first derived in [30, 1] for constant-curvature 
spaces SO(N + l)/SO(N) ~ S N and subsequently generalized in [3] to the spaces 
G/SO(N), covering all examples in which the bi-Hamiltonian structure is a <3(iV — 
l)-invariant (namely, HZ = 0(N - 1) C SO(N)). The proof of Theorem 4.5 in the 
general case is similar to the special case proven in [30] where t)\\ — so(N — 1) ~ 
r n-i aR a-i^ ^ _ r n-i such that [fj^^n] c f) ± and [Ju,fn] = f)||. Note these 
Lie bracket relations correspond to the geometrical properties (39) of a Riemannian 
parallel frame [3]. There are two main differences in the general case. Firstly, the 
operators 7i, J have an additional term arising from the Lie bracket structure 
[f)j_, f)j_] n t)± ^ as permitted by the more general algebraic notion of a parallel 
frame in Definition 3.3 and Theorem 3.5. Secondly, these operators are formulated 
here in a manifestly iJ|*-invariant fashion, using just the intrinsic Lie bracket on 
i) = f)|| © f)j_, as opposed to a more cumbersome R JV_1 -component form used in [30]. 

4.3. Proof of cosymplectic property for H. A linear operator H that acts on 
()x-valucd functions in the x-jet space J°° of the flow variable it is a Hamilton- 
ian cosymplectic operator [29, 15] iff its associated Poison bracket (80) is skew- 
symmetric and obeys the Jacobi identity. 

Skew-symmetry (81) is equivalent to H being skew-adjoint with respect to the 
natural inner product J (/i, f2)dx on fjj_-valued functions / on J°° induced by 
the Cartan-Killing inner product on hj_. This condition can be formulated in a 
pointwise manner on J°° as follows: H is skew-adjoint iff 

(f,H(f)) = modulo a total x-derivative (87) 

(namely, (f,H(f)) = D x T for some scalar function T). 

Proposition 4.7: The linear operator H = D x + [u,-]± — [u,D~ [u, -]||] in 
Theorem 4-5 is skew-adjoint. 

Proof. The l.h.s. of (87) is given by 

(f,DJ + [ U J] ± -[ U ,D- I {uj] ll ]} = (f^uJ^-if^^D-'iuJ}^) (88) 

since (f,D x f) = —\D x \f\ 2 . The remaining two terms can be simplified by means 
of the identities 

(A, [B,C±]) = cyclic, (A, [B,C\\]) = (C h [A,B]), (89) 
for functions A, B with values in f)j_ and C in (). Thus 



(/, [u,f]) = -(u, [/,/])=0 



(90) 
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and 

(/, [u,A(f)}) = -(A(f), [«,/]„) = \D x \A(f)\ 2 = (91) 

where A(f) = D~ 1 [u, f] \\ . Hence the l.h.s. of (87) reduces to a total x-derivative, 
which completes the proof. 

For a skew-adjoint operator TL 1 the Jacobi property (82) is well-known to be 
equivalent to the vanishing of the Schouten bracket of 7i with itself [15]. The 
simplest approach to verifying if an operator has vanishing Schouten bracket is the 
calculus of multi-vectors developed in [29] which will be adapted to the ()x-valued 
setting here. 

First, on x-jet space J°°, introduce a f)j_ -valued vertical uni-vector 9 which is 
dual to du regarded as a f)_L-valued 1-form, in the sense that —(6,du) = 1 where 
the pairing denotes the usual contraction between vectors and 1-forms, combined 
with their Cartan-Killing inner product on t). Then x-dcrivatives of 9 are introduced 

2 . 2 

similarly as duals to du x = D x du, du xx = D x du, etc. via 9 X — D x 9, 9 XX — D x 9, etc. 
such that —(9k x , dui x ) = Ski (with k, I = 0, 1, 2, . . .). Formal sums of antisymmetric 
tensor products of these vertical uni-vectors will define vertical multi-vectors, for 
instance 9A9 and 9/\9/\9 with respective values in Aj? ± and Ajj^ . More particularly, 
a combined wedge product and Lie bracket of 9 with itself yields a [)-valued vertical 
bi- vector 

[M = [M]± + [M|| (92) 

while the vertical tri-vector 

[[e,6],0]=0 (93) 

vanishes due to the Jacobi identity on f). 

It is useful to let the basic uni-vector 9 have the formal action of a vector field 
on f)^-valued functions f(x, u, u x , . . .) on J°°, 

Sof^^ekxdf/duk*. (94) 

0<fe 

Likewise, let 9 act on differential operators T> = X)o<fe by a formal Lie deriv- 

ative 

S V := 5^(J e o fc )D* (95) 

0<fe 

in terms of coefficients ak(x,u,u x , . . .) given by functions on J°° . In particular, 
SeD x = 0. Finally, the action of 9 can be naturally extended first to vertical 
uni-vectors 9 via defining 

5 e 9 := (96) 

so 5e9k x — 0, and then to vertical multi-vectors by letting 5$ act as a derivation. 

In this formalism the Schouten bracket of TL with itself vanishes if and only 
if J (6 A dx = where (• A •) denotes the Cartan-Killing inner product 

combined with the wedge product of f)j_-valued 1-forms on J°°. Note, in particular, 
(9 A 9) = 0. An equivalent pointwise formulation is that (9 A 8-n(g)H(8)) = D X T 
(for some scalar function T) or simply 

(9 A 5t-i(q)H{9)) = modulo a total x-derivative. (97) 

This condition is equivalent to the Jacobi identity (82) holding for the Poisson 
bracket {•,•}«• 
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Proposition 4.8: The linear operator H = D x + [u,-]± — [u, D~ [u, -]||] in 
Theorem 4-5 has vanishing Schouten bracket and hence is cosymplectic. 
The following multi-vector identities will be used in the proof: 

[A(9),B(9)] = [B(9),A(9)}, (98) 

[A(9) : [B(9),C{9)]] + cyclic = 0, (99) 

and 

{A(9)AB(9))=-(B(6)AA(6)), (100) 

(A(9) A [B{6),C{9)\) = cyclic, (101) 

for linear operators A,B,C from fjj_ into h; also, 

(A(9) A [u, •]) = —(u, [A{6), •]) = —([u, A{6)] A •), (102) 
[[«, + [[«, 4(0)], = [«, (103) 

and in particular 

[[u,A(fl)],A((9)]] = l[«,[A((9),A(<?)]]. (104) 

Proof. Consider 

5nmH(6) = [H(9),9] ± - [H{9),D x \u,6]\ti - [u, D~ X \H{6), 6]\\] (105) 

which expands out to a ^-valued cubic polynomial in u. There is a single inho- 
mogeneous term, [D x 9,9]_\_. This yields in the l.h.s of (97): 

(0A[D X 6,O]) = ±D X {9A[6,9])=O (106) 

by use of (101). 

Next the linear terms in (105) are given by 

- [DJ, D-^u, 9} || ] - [«, D-\D x e, 0\\\] + [[u,0\±,0\±. (107) 
In the l.h.s of (97) the first term of (107) yields: 

-{9A[D x 9,D- 1 [u 7 9] ll }) = —(D-'^e^AD^e}) = |<[^,^]|| A [9,9]) (108) 

from (101) and (98). Similarly the middle term of (107) simplifies: 

-(9 A [u, D-^DJ, 9} ||]) = («, [9, D-^DJ, 9} „]) = - [0, [0, 0] J) 

= -^<[M]A[M±) (109) 
from (102), (92) and (93). The last term of (107) directly gives: 

(9 A[[u,e} ± ,9}) = ([u,6} ± A[9,9}). (110) 
Hence the linear terms in the l.h.s of (97) combine into 

= \{[u, 9] || A [9, 9}) + \ <[M] ± A [9, 9}) = ±([u,9]A [9, 9]) = \ (u, [9, [9, 9]]) =0(111) 

by (102) and (93). 

The quadratic terms in (105) are given by 

-[MU^Mh] - [[^D-^e^eu [^^[MU,^]. (112) 
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Write A{8) = Z?^ 1 [u, 0] y and combine the first and second terms in (112) by the 
Jacobi relation (99) to get [[A(8),6],u]±. This yields in the l.h.s of (97): 

(8A[[A(8),8],u]) = (A(8)A[[u,8],9}) (113) 

from (102) and (101). For the third term in (112), we use (102) to simplify the l.h.s 
of (97): 

-<6»A[«,D; 1 [[«,e]x,fl]||]) - (MIAZ^IMU^h) = -{A(8)A[[u,9U,8}). (114) 
This term now combines with the previous one, giving 

ee (A(8) A [[«,%*]) = (A(8) A [D x A(0),6]) = \{9 A D x [A(9), A(8)} ± ) (115) 

from (101). But the Lie bracket relation (66) implies [A(8), A(8)]j_ vanishes and 
hence so does the combined quadratic term (115) in the l.h.s of (97). 
Finally, there are two cubic terms in (105): 

[[u,A(8)],A(9)} + [u,D x 1 [[u,A(0)},9] ]l ] (116) 

where A{8) = D~ [u, 8]\\ again. The second term of (116) simplifies in the l.h.s of 
(97): 

(8A[u,D- 1 [[u,A(8)],8] ll })^-([u,8}AD- 1 [[u,A(8)},8} ll ) 

ee (A(8) A [[«, A(8)],8]) = (8 A [A(9), [«, A{8)]]) (117) 

by (102) and (101). This term combines with the first term of (116) in the l.h.s of 
(97): 

2(0 A [[u,A(9)},A(8)}) = (9 A [u, [A(8), A(9)]\{\) = -(D X A(8) A [A(9),A(9)}) 

= -±D x (A(9)A[A(9),A(9)})=0 (118) 

by means of (104), (101) and (102), in addition to [A(9),A(9)} ± = 0. 

Thus the l.h.s of (97) vanishes modulo total x-derivatives, which completes the 
proof. 

4.4. Proof of symplectic property for J. A linear operator J that acts on de- 
valued functions in the x-jet space J°° of the flow variable u is a symplectic operator 
[15] iff the associated bilinear form (84) is skew and closed, namely u>(X,X) = 
and doj(X, Y, Z) = 8 x u{Y, Z) + cyclic = for all vector fields X, Y, Z on J°°. 

The conditions on w have a straightforward pointwise formulation on J°° in 
terms of 

wj(h±,h£):=-(h±,J{h£)) (119) 
where h 1 - denotes a f)^-valued function. It will be convenient to work equiv- 
alently with m^-valued functions h± = ad(e a ,) _1 /i- L along with the operator 
JC\ m ± = -ad(e x ) l 7ad(e x ) as follows: 

u>j(hi,h£) - (ad(e x )/ lJ . 1 ,ad(eJ- 1 /C(/ lJ . 2 )) = -(h ±1 ,IC(h ±2 )). (120) 

Note, here, ad(e x ) is skew-adjoint with respect to the Cartan-Killing inner product 
on hj_ ® mi as a consequence of the elementary identity 

(ad(e x )A, B) = (e x , [A, B]) = -(A, ad(e x )B) (121) 

for functions A, B with values in \j± © mj_. 
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To proceed, skew-symmetry of u> will hold iff uij is skew, which is clearly equiv- 
alent to JC being a skew-adjoint operator: 

(/, JC{f)) = modulo a total rr-derivative (122) 

for all tn_i_-valued functions / on J°°. This condition is formally identical to the 
skew- adjoint property established earlier for the operator H = /C|j, x in Proposi- 
tion 4.7, where, recall, JC is the operator given in (74). As a result, the isomorphism 
mj_ — f)± provided by ad(e a; ) implies that the steps in the proof of Proposition 4.7 
carry through with the function / now being m_i_-valued instead of f)j_ -valued. 

Proposition 4.9: The linear operator J — — ad(e x )~ 1 (D x + [u,-]± — 
[u, D~ X [u, -]||]) ad(e 2 ,) _1 in Theorem 4-5 is skew-adjoint and hence the bilinear form 
= J u>j(-, -)dx is skew- symmetric. 

Next, for a skew-adjoint operator J , the closure property du} = is equivalent 
to having uij satisfy S h ±cjj(h 2 ,h 3 L ) + cyclic = D x T for some scalar function T. 
This condition will hold iff: 

(h ±3 ,5 hi fC(h ±2 ))+ cyclic = (123) 

for all m_i_ -valued functions h±i, with hf- = &d(e x )h±i (i = 1, 2, 3). 
Hereafter it will be useful to let 

(-,.):= [ad( e j,-] (124) 

which satisfies the following identities: 

(A,(B,C))=-(C,(B,A)) (125) 

and also 

(A,B)||=(B,A)||, (126) 

(A, B) ± - (B, A) ± = ad(eJ[A B], (127) 

(A, (B, C)±) + cyclic = 0, (128) 

for functions A, B, C with values in mj_. Note the latter two identities are readily 
established from the property that ad(e a; ) acts as a derivation on [•, •] in and 
annihilates my . 

Proposition 4.10: The skew bilinear form = J u)j(-,-)dx is closed and 

hence the linear operator J = —a,d(e x )~ 1 (D x + [u, -]± — [u, D" 1 ^, •] ||])ad(e a ,) — 1 in 
Theorem 4-5 is symplectic. 

Proof. Consider, from (74) and (124), 

= (h±i,h±2}± - (h±i,D~ [u, h± 2 }\\) - [u,L>^ 1 (/i_li,/i_l 2 )||] 

(129) 

which is a mj_ -valued linear polynomial in u. The inhomogeneous term yields in 
the l.h.s of (123): 

(h ±3 ,(h ±1 ,h ±2 )±) + cyclic = (130) 

from (128). Write A(h±) = D^ 1 [u, h±] || . Then in the l.h.s of (123) the middle term 
in (129) is given by: 

-{h±3,(h±i,A(h± 2 ))) + cyclic = (A(h ±2 ),(h ±1 ,h ±3 )\\) + cyclic (131) 
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from (125). For the last term in (129), the l.h.s of (123) simplifies: 

-(h±3, [u, D~ (h±i,h± 2 )\\}) + cyclic = ([u,h± 3 ],D~ (h±i,h± 2 )\\) + cyclic 
= -(A(h ±3 ), {h ±1 ,h ±2 )) + cyclic = -(A(h ±2 ), (h ±3 , h ±1 )) + cyclic (132) 

by cyclic rearrangement. Note A(hj_)± = 0. Hence, from (126), the terms (131) 
and (132) cancel, 

= (A(h ±2 ), {h ±1 ,h±s)\\ - (h ±3 , h ±1 )\\) = (133) 
which completes the proof. 

4.5. Proof of compatibility of 7i,J. Hamiltonian cosymplectic and symplectic 
operators TL, J are said to be compatible (i.e. comprise a Hamiltonian pair) if every 
linear combination c\H + c 2 J~ 1 is a Hamiltonian operator. Equivalently [29, 15], 
compatibility holds whenever Ti and J^ 1 have a vanishing Schouten bracket. 

The latter condition can be formulated using the same calculus of vertical multi- 
vectors as employed in the proof of Proposition 4.8: 

(8 A {5 H {e)J^ 1 {G) + Sj-i(0)H(6))} = modulo a total x-derivative (134) 

where 8 denotes the vertical uni-vector dual to du in the x-jet space J°° of u. It 
is possible to simplify the l.h.s. of (134) by noting 5n(e)J^ 1 = —J~ 1 {o~h{b)J)J~ 1 
and then replacing 8 by J{8) followed by using the skew-adjoint property of J to 
get 

(0 A S n( jm J(8) + J{8) A (5 e H)J(6)) = 0. (135) 

Now, as the uni-vector defined by J{8) obviously enjoys the same formal properties 
as 8 itself, this establishes that (135) is equivalent to condition (134). 

Lemma 4.11: The Schouten bracket ofTL.J vanishes iff (135) holds. 

For verifying if H,J satisfy (135), the earlier multi-vector identities (98)-(102) 
will be useful, in addition to the following identities: 

(A(8), (B(6), C{0))±) + cyclic - 0, (136) 

(A(8),B(8)) ll =-(B(8),A(8)) h (137) 

and in particular 

(A(8) 7 A(8)) ll =0 7 (138) 

for linear operators A, B, C from t)± into f). Here, (•, •) = [ad(e x )- ; •] as before. 

Theorem 4.12: The cosymplectic and symplectic operators Ji^J in Theo- 
rem 4-5 have vanishing Schouten bracket and hence comprise a Hamiltonian pair. 

Proof. It will be useful to introduce the uni- vectors 8 := a,d(e x )~ 1 8 and 8 := 

ad(e a ,)" 1 0, with values in and f)j_, having the same formal properties as the 
(^-valued uni-vector 8. (In particular, 8 = &d(e x )8 — a,d(e x ) 2 8.) To begin, note 

(8 A b~H(j(8))J{8)) = (8 A o~u(j{6))fc-{8)) due to ad(e 2; ) being skew-adjoint with 
respect to the Cartan-Killing inner product. Now consider 

WwW) = -wm^wM - [u,D?[nJV))Aw] + wmwh 

(139) 

and write A{8) := D x 1 [u, 8} \\ . In the l.h.s. of (135), the first term in (139) yields 



-(6*[H(J(0)),A(e)]) = 



-(A(8)A[8,H(J(8))}) 



(140) 
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and the second term in (139) simplifies: 

-(6A[u,D?[H(J(e)),§]\\]) = -(A(6) A [6,H(J(6))]). (141) 

Hence these two terms add to give, from (73), 

-2(A(9) A [6, D x J{6)]) - 2(A(9) A [[«, J(6)] ± ,6]) + W~9) A [[u, B{6)],9]) (142) 

where B(6) := D' 1 ^, J{9)}\\. Next, in the l.h.s. of (135), the third term in (139) 
expands out to yield 

(§A[H{j(0)),§]) = ([9,9]ad x j(6))+([9,9]a[u,j(6)] ± )-([9,6]a[u,b(9)]). (143) 

Note the six terms given by (142) and (143) are linear in J (6). Six similar terms 
come from (J(6) A (5eTL)J(9)) as follows. First consider 

(5 $ H)J(9) = -[6,D- 1 [u,J(6)] ll ] - [u,D- 1 [9,J{9)\] + [9,J{9)] ± . (144) 

In the l.h.s. of (135), the second term in (144) simplifies via (101) 

-(J(9) A [u 7 D- 1 [9,J(9)} ll ]) eee -(B(9) A [6, J (9)]) (145) 

which then adds with the first term in (144), yielding by (137) and (73), 

-2(B(9) A [e,J(6)}) = -2(B(9) A (0,J{0))\\) = -2{B{9) A [JC(0),0]). (146) 

This expands out to give three terms: first 

-2{B{9) A [Dj,0]) ee ([u,J{9)]\\ A [W]) (147) 

by (98); next 

-2(B(9) A [[u,0]_l,0]) = -2(B{9) A [M],#]||) = ~{B{0) A [u,[O,0\]) (148) 
by [[u, 0]\\, 0]\\ = from the Lie bracket relations (67) followed by (104); third 

2(B(9) A [[u, A(9)],9}). (149) 
The last term remaining in (144) can be rearranged in the l.h.s. of (135): 
(J (9) A [9,J(0)}) - (J(9) A (0,J(0))±) 

= -(9 A (J (9), J (9)) ± ) - (J(9) A (J(9),9)±) - 2{J(9) A [9,K(§)])(150) 
via (136) and (101). This term expands out to give three terms: 

(J(9) A Dj, §]) - 2(J(9) A [9, [u, A{9)]]) + 2(J(9) A [9, [u, §\ ± ]) (151) 
via (98). 

Now the term (148) cancels the last term in (143) by (102), while the first terms 
in (151) and (143) reduce to DJJ(9) A [6,6]) = 0. 

The remaining terms in (151) and (143), along with the three terms in (142) and 
the two terms given by (147) and (149), constitute a cubic polynomial in u with 
coefficients linear in J (9). The linear u-terms in this polynomial consist of: 

2{J(9) A [6, [u, §] ± ]) = (J(9) A [u, [9, 9]]) - 2{D X A{9) A [6, J{9)]) (152) 

via (104) and (101); plus 

([6,6] A [u, J (9)]) (153) 

and 

-2(A(6)A[9,D x J(6)]). (154) 
When combined, these four terms reduce to 

= 2{A{0)A[Dj,J{0)]). (155) 
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Next, the quadratic u-terms in the polynomial yield 

-2(J(0)A[0,[u,A(0)}])-2(A(0)A[0,[u,J(0)] ± ]) = -2(0A[u,[A(0),J(0)}]) (156) 

via (101) and (103), with [0, [u, J(0)]\\]\\ = by the Lie bracket relation (67). 
Similarly, the cubic terms combine to produce 

2([[u,A(0)},0] A B{0)) + 2(A(0) A [[u,B(0)],0]) = 2(0 A [u, [A{0),B(0)]) 
= -2(D X A(0) A [A(0),B(0)}) ee ([A(0),A(0)} A [«, J{9)]) (157) 

where the middle step uses [A(0), B(9)]± = due to the Lie bracket relations (66). 
To proceed, we make use of the identity 

J(0) = ~ad(eJ- 1 /C(0) = -D x C{0) + [u, A(0)} (158) 

with C(9) := ad(e a ,)~ 1 [u, 0]± and u :— ad(e a ,) _1 u where the last term in (158) 
arises from [u, A(0)\ — &&(e x )[u, A(9)\ as ad(e a; ) annihilates A(9). The terms (155), 
(156), (157) expand out to become a quintic polynomial in u as follows. There is a 
single linear term, coming from (155): 

-2(A{0) A [Dj,D x 9]} = -2(A(9) A (D X 6,D X 0)\\) = (159) 

by (138). Next, one quadratic term comes from (155): 

-2(A(9) A [Dj,C(9)}) = -2(A(9) A (D X 6,C0))\\) = 2(A(0) A [[u,0] ± ,D x 0]) 

= 2(D x 0A[A(0),[u,0]}) (160) 

via (138) and (101) combined with [A(9), [u, 0]\\]± — by the Lie bracket relations 
(66). This term (160) now cancels the quadratic term given by (156), 

2(0 A [u,[A(O),D x 0]}) = -2([u,0] A [A(9),DJ}} (161) 

through (101). We then have only cubic and higher degree terms remaining to 
consider. One cubic term comes from (156): 

2(0 A [u, [A(0),C(0)}}) = -2(A(0) A [\u,0],C(0)}) - -2(A(0) A {C(0),C(0))\\) = 

( 162 ) 

by (138), with [[u, 0]\\, C(9)]^ = by the Lie bracket relations (67). Another cubic 
term arises from (155): 

2(A(0) A [Dj, \u,A(0)}]) = (D x 0A[u, [A(0), A{0)]]) = ([A(0),A{0)\ A (D x 6,u) n ) 

(163) 

via (101) and (104). A similar cubic term is given by (157): 

([A(0),A(0)} A [u,D x 0]) = -([A(9),A(0)] A (u,D x 6) n ). (164) 

Then (163) and (164) cancel by (126). There is a similar cancellation of quartic 
terms coming from (156) and (157): 

-2(0A[u, [A(9),[u,A(0)}}]) = ([A(6), A(0)]A[[u, 0] ± , u}) = ([A(0),A(0)]A{C(0),u)\\) 

(165) 

and 

-([A(0),A(0)} A [u,C(0)}) = -([A(0),A(0)] A (u,C(0)) l{ ). (166) 
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As a result, the only term remaining is quintic, which comes from (157): 

([A(e),A(§)}A[u,[u,A(e)}}). (167) 
We now employ the Jacobi identity combined with the identity 

[u,[«,A($)]]|| = («,[«, = ([fi,A(e)],«)|| - [[u,A(0)],fi]|| (168) 

to get 

[ U ,[«,A(fl)]]|| = ![[«,«], (169) 

Hence (167) simplifies to 

i([A(0),^)] A [[«,«], - ±<[u,fi] A [[^),A(0)],A(0)]> = (170) 

by (99). Consequently, all terms in the l.h.s. of (135) have been shown to vanish 
modulo a total rr-dcrivative, thereby completing the proof. 

5. Bi-Hamiltonian hierarchies of non-stretching curve flows and 

soliton equations 

In any Riemannian symmetric space M — G/H, the frame structure equations 
of non-stretching curve flows ~/(t,x) geometrically encode a group-invariant bi- 
Hamiltonian structure for the induced flow on the components of the principal 
normal vector along the curve. This structure looks simplest as stated in Theo- 
rem 4.5 using an if -parallel moving frame formulation (cf Definition 3.3) expressed 
in terms of a m- valued linear coframe e and its associated ^-valued linear connection 
oj, obeying properties (i) to (iv) summarized in Lemma 4.1. Such a framing of j(t, x) 
arises directly from the g- valued zero-curvature Cartan connection s w = e+uo in the 
Klein geometry associated with G viewed as an .ff-bundle over M , where g = m© f) 
is the corresponding Lie algebra decomposition of G and m ~g T x M. In particular, 
e and u> provide a soldering of the Klein geometry of G/H onto the Riemannian 
geometry of M, such that the frame structure equations coincide with the Car- 
tan equations for torsion and curvature (cf Theorem 2.3 and Corollary 2.4) of the 
Riemannian connection V t , V x on the two-dimensional surface of the flow in M. 

The bi-Hamiltonian structure consists of compatible Hamiltonian (cosymplcctic 
and symplectic) operators Ti and J, which yield a (hereditary) recursion operator 

Tl = HJ = -(/Cad(e x )- 1 ) 2 

= -(^(eJ-'D^iu^dieJ-'-h-^D-^u^dieJ- 1 -]^}) 2 (171) 

with respect to the flow variable u. Due to the ii-parallel property of e and oj, 
e x = &\lx is a fixed unit vector in m representing the moving frame components 
of the tangent vector ^ x along the curve, while u — oj x = oj\j x is related to the 
moving frame components of the principal normal vector V x j x along the curve 
through the corresponding m^-valued expression [u,e x ] — — ad(e x )u = e\V x ^ x . 
Note ad(e a ,) _1 gives an isomorphism between the perp spaces f)j_ and of the 
centralizer subspaces f)|| and rrt|| with respect to e x in m © t) = g, with ad(e x )f)|| = 
ad(ea;)m|| = 0. Hence the linear operator /Cad(e a .) _1 maps f)j_-valued functions into 
mj_-valued functions, and vice versa. Thus the adjoint of the recursion operator is 
given by 

TV - 1 7W = -(ad(eJ- 1 /C) 2 

= -(adtej- 1 ^ + ad(ej- 1 [« ) • ]u - ad^)" 1 ^,^ 1 ^, • ],|]) 2 . (172) 
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On the x-jet space J°° = (x, u, u x , u xx , ■ ■ ■) of the flow variable u, the x- 
translation vector field d/dx is an obvious symmetry of the operators TL,J and 
hence of the recursion operator 1Z. As a consequence, from general results due to 
Magri [23, 24] on Hamiltonian recursion operators [29], 1Z will generate a hierarchy 
of commuting Hamiltonian vector fields with respect to the Poisson bracket, where 
the hierarchy starts with the evolutionary form of the vector field d/dx. More- 
over, an associated hierarchy of involutive covector fields arises from the canonical 
pairing provided by the symplectic 2-form. 

Theorem 5.1: A commuting hierarchy of Hamiltonian vector fields ■ d/du 
is given by the t)±-valued functions 

hf n) :=TZ n (u x ), (n = 0,1,2,...) (173) 

satisfying the Poisson bracket relation 

6 hHn) - = {;rt n) } H (174) 

for some Hamiltonian functionals S)^ = J H^ n \x,u,u x , . . .), with respect to the 
Hamiltonian operator H. The variational derivative of each such functional yields 
a dual vector field • du given by the f) ±-valued functions 

wf n) := SH {n) /8u = K* n (u), (n = 0, 1, 2, . . .) (175) 

forming a hierarchy of involutive variational covector fields. In particular, at the 
bottom of these hierarchies, H^ — —^(u,u) is the Hamiltonian for the vector field 
given by /itL = u x corresponding to the generator of x -translations in evolutionary 
form, whose dual covector field is given by = u. 

Corollary 5.2: These hierarchies are related by the Hamiltonian operators 

^f n+ l)=^f n) ), Kn)= U Hn)) ( 176 ) 

with {'OJ'hy h(~ n )) = modulo a total x-derivative, since Ti., J are skew-adjoint. 

The form of these operators under a scaling of the variables x and u shows that 
both hierarchies possess the mKdV scaling symmetry x — ► Ax, u — ► \~ 1 u. Thereby 
hh n ^ and H^ each have scaling weight 2 + 2n, and hence has scaling weight 

l+2n. From variational scaling methods in [2], it follows that H^ can be expressed 
explicitly in terms of h^ or \ : 

(l + 2n)If (n) =-D- 1 {u,hf n) ) = -D-'^H^tn)))- ( 177 ) 

The derivation of this formula is fairly simple. Let S s u = —xu x — u = —D x (xu) 
and 8 S H = —xD x H — pH be the evolutionary form of the mKdV scaling symmetry 
generator on u and H(x,u,u x , . . .), with respective scaling weights 1 and p ^ 1. 
Then we have 

6 s H=(l-p)H (178) 
modulo a total x-derivative. On the other hand, 

5 S H = -{8 s u,8H/5u) = (D x (xu),w ± ) = -x{u,D x vj^) = -x(u, H{vj^)), (179) 

where the last line uses the identity (u, [u, Afa- 1 )]) = (A(w A -) 1 [u, u]) = holding 
for any linear operator A. Finally, the identity (u, 7i(zo )) = D x D x (ujHinj- 1 )) 
yields 5 S H = D x 1 (u, TL{w^)) modulo a total x-derivative, from which we obtain 
H=(l-p)^D- 1 (u,H(^)). 
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5.1. Group-invariant bi-Hamiltonian soliton equations. The entire hierar- 
chy of vector fields (173) has two compatible Hamiltonian structures Tt and J~ x : 

hf n) =H(SH {n) /Su) = J- l (5H (n+1) /Su), (n = 0, 1, 2, . . .). (180) 

The operator £ := 1ZH provides an alternate Hamiltonian structure 

hf n) = S{6H {n ~ 1) /Su) (forn>0), (181) 

with Ti and £ comprising an explicit Hamiltonian pair (of cosymplectic operators). 
As a main result, this bi-Hamiltonian structure (180) produces a hierarchy of bi- 
Hamiltonian flows on u as given by linear combinations h 1 - = h^ n+1 ^ + of the 
vector fields (173) according to the form of the flow equation (72). 

Theorem 5.3: The flow equation (72) on u yields a hierarchy of bi-Hamiltonian 
evolution equations 

u t =H{SH (n > 1) /Su) = J-^SH^'V/Su) (182) 

(hereafter called the +l + n flow) for n = 0, 1, 2, . . ., where i^™' 1 ) : = +H^ n ~ 1 " 1 . 
Each of these t)±-valued evolution equations is invariant under the equivalence group 
Ht C H* consisting of rigid (x -independent) linear transformations Ad(ft, _1 ) that 
preserve e x . 

These evolution equations (182) concretely describe group-invariant multicom- 
ponent soliton equations when u is expressed by real- valued components in an ex- 
pansion with respect to any fixed basis (11) and (22) for the vector space t)± ~ mj_. 
These components will naturally divide into irreducible representations u a un- 
der the action of the group ify, corresponding to a direct sum decomposition 
f)_L = ® Q f)" where each subspace ()" is an eigenspace of Ad(if|*). The specific 
algebraic structure of these eigenspaces (which will depend on the structure of 
the group _ff|) defines the algebraic content of the variables u a : e.g. scalars, real 
vectors, complex vectors or vector pairs, scalar- vector pairs, etc.. For instance, 
in the example M = SO(N + l)/SO(N) studied in [3], the Lie algebra of i/j is 
t)|| = so(N — 1) c f) = so(N), whence \)± ~ ]R W_1 ~ with this vector space 
being irreducible (i.e. a = 1). Thus, in this case, u G R^ -1 is algebraically 
identified as a single real vector variable, and the hierarchy of vector evolution 
equations (182) describes SO(N — l)-invariant multicomponent soliton equations 
for u = (ti 1 , . . . , u w_1 ). 

In general the bi-Hamiltonian evolution equations (182) can be rewritten in terms 
of mj_ -valued variables more geometrically related to the flow ~/(t,x) as follows: 

N = V xlx ^v:= -ad(e x )u (183) 

represents the principal normal vector along the curve; 

(7t)± ^ h± := ad(e a; )" 1 /i- L , {lt)\\ ^> h\\ (184) 

represent the algebraically perpendicular and parallel parts of the flow vector; 

V t 7a; <-» zu± := — ad(e a; )n7 ± = — a,d(e x )(zu ± + ro") (185) 

represents the principal normal vector relative to the flow. These geometrical vari- 
ables inherit a natural Hamiltonian cosymplectic and symplectic structure from the 
results in Theorem 5.1 and Corollary 5.2 expressed with respect to the flow variable 

v := &d(e x )~ 1 u = X^ x {v), X = -&&{e x f (186) 
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in the following way. To begin, the induced flow equation takes the form 

v t = H(w JL ) + h ± , w JL = J(h 1 _) (187) 

where 

H = -a,d(e x )- 1 )C\ t , ± a,d(e x )-\ J = K\ m± (188) 

are compatible Hamiltonian cosymplectic and symplectic operators given in terms 
of 

K = D x + [ad{e x )v, ■ ] ± - [ad{e x )v, £>; 1 [ad(e x )u, ■ ]y] ± (189) 

from (74) and (186). 

Proposition 5.4: On the x-jet space of geometrical flows with respect to v, 
h± ■ d/dv is a Hamiltonian vector field 

S h ± (•) = {-,?)} fl = - J (6(-)/Sv, H(6f)/6v))dx (190) 

with respect to the Poisson bracket {■, -}^, and zu± ■ dv is a variational covector 
field 

w± = 5$/5v (191) 

satisfying the duality relation 

(n7 ± • dv)\ ■ =u(- ,h±- d/dv) j = — J (■ , J(h±))dx (192) 

with respect to the symplectic 2-form uj(-, ■) j. 

Theorem 5.5: There is a bi-Hamiltonian hierarchy of evolution equations on 
the principal normal components (183) in a H -parallel moving frame, represented 
by 

v t =H(5H {n < 1) /5v) = j- 1 (8H {n+1 - 1) /8v), (n = 0, 1,2, . . .) (193) 

with Hamiltonians H^ n ^ ■= + H^^ 1 ^ expressed in terms of the mj_-valued 
variables (v,v x , . . .) via (186) and (177). Correspondingly, the normal flow compo- 
nents 

h™=it n (v x )=H(m™), (n = 0,1,2,...) (194) 
and the principal normal components with respect to the normal flow 

w ( [ l) =1Z* n (X(v))=J(h^- 1) ), (n=l,2,...) (195) 

represent hierarchies of commuting Hamiltonian vector fields ■ d/dv and invo- 

the recv 

TZ = HJ = -(ad(e a; )- 1 /C) 2 (196) 



lutive variational covector fields td^ ■ dv generated by the recursion operator 



and its adjoint 

TZ* =JH = -(lCad(e x )- 1 ) 2 . (197) 

This version of the flow equation and its bi-Hamiltonian structure is closest to 
the notation used in [3, 4]. Note the invariance group of the operators H and J is 
i?| C H* ~ Ad(H), and the scaling weights of h^ and are 2 + 2n and 1 + 2n 
under the mKdV scaling group x — > Xx, v — > A _1 w. 
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Through the geometric correspondences (183)— (185), the evolution equations 

(193) produce a bi-Hamiltonian hierarchy of non-stretching curve flows in M — 
G/H. The Hamiltonians are related to the tangential part of the flow by 

fl(7t,7x) - -(h (n \e x ) = {l + 2n)H {n) (198) 
arising from (70) and (177), where 

/I,, = D-^vMieMw, H {n) = --l_ D - l {v,X{h^)). (199) 
Remarks: 

(1) The infinitesimal invariance group of the hierarchy (193) can be identified 
with the Lie subalgebra f)|| C f) given by the centralizer of e x in the Lie algebra of 
the group H. 

(2) The bi-Hamiltonian structure (176) of the hierarchy depends on the constant 
unit vector e x in m ~ g/f) determined by the framing of the curve flow j(t, x). 

(3) Two curve flows encode equivalent (isomorphic) bi-Hamiltonian structures iff 
e x lies on the same orbit of the linear isotropy group H* in m for both flows, where 
all such orbits are parametrized by the unit-norm elements in the Weyl chamber 
a*(m) (cf Proposition 3.2) contained in a fixed maximal abelian subspace a C m. 

(4) Without loss of generality, in a given equivalence class, a representative e x 
can be chosen such that the linear map ad^) 2 is diagonal on m^. 

Consequently, every Riemannian symmetric space M = G/H will possess dis- 
tinct bi-Hamiltonian hierarchies of non-stretching curve flows characterized by the 
distinct elements in a canonical Weyl chamber o*(m) C m ~ g/f) ~ G T X M . For 
the example M = SO(N + 1)/ SO(N), it is well-known that every maximal abelian 
subspace a C m = so(N+l)/so(N) ~ R N has dim a = 1, whence the Weyl chamber 
a*(tn) has a single unit-norm element. As a result, up to equivalence, there is a 
single bi-Hamiltonian hierarchy of curve flows in M = SO(N + 1)/ SO(N). 

The equations of the curve flows on j(t, x) in a given hierarchy in M = G/H can 
be readily derived in an explicit form from the geometric correspondences (183)- 
(185) by means of Proposition 4.4. 

5.2. mKdV flows. The +1 flow on v will now be derived in detail. This will 
require working out the first higher-order Hamiltonian structures in the hierarchy 

(194) and (195). At the bottom of the hierarchy, note 

hf=v x , (200) 
ff (0) = -{h^\e x ) = -D-^v^ihf)) = -l|ad(e>| 2 , (201) 
wf =SH (a) /Sv = X(v), (202) 

along with 

h\V = D-\vM{e x )h^ ],| = -i[ad(e>,«]||, (203) 
where, recall, X = —a,d(e x ) 2 . The next order in the hierarchy is given by 
SH {1) /Sv = vu ( l ] = J(hf>) = TZ*(X(v)) 

= v xx + [ad(e x )v,v x -^[ad(e x )v,v]\i], (204) 
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h { l ] = H{w ( l ) )=n{v x ) 

= ^ _1 Kxi) + X~ 1 ([ad(e x )v,v x - ^[ad(e x )v, v]^]±) x 



-ad(e x ) 1 [ad(e x )v, ad(e x ) 1 (v xx + [ad(e x )v, v x ]±) 



[ad(e x )u,u]||]] - [v, [v,v x + ^[ad(e x )v,v]±\u], (205) 



= [aA(e x )~ 1 v xx ,v]\i - ^[a.d(e x )~ 1 v x ,v x }\\ + [ad(e x y 1 v, [ad(e x )v, v x ]]\\ 
+ [v, [ad(e x )- 1 v x ,ad(e x )v]]\\ + ^[v, [v,3[ad(e x )v,v]\\ + [ad(e x )u,u]j.]]|| 
-^[&d(e x y 1 [ad(e x )v,v} ± , [ad(e x )v,v] ± \ l \, (206) 

H (1) = -(hf\e x ) = -D- 1 (v,X(h^)) 

= ~\v x \ 2 + ^(v x ,[ad{e x )v,v] ± ) + ^\[ad{e x )v,v} ]l \ 2 . (207) 

The derivation of (204) is simple: The local part of J on (200) yields 

D x hf + [ad(e x )v, h^U = v xx + [ad(e x )v, v x } ± , (208) 
while the nonlocal part of J applied to (200) reduces to 

[ad(e x )v, /ij 0) ] = ~[ad(e x )v, [ad(e x )v,v]\\]- (209) 

For the derivation of (205), the following identity will be needed: 

[A, (B, C)] || - [C, (B, A)] || = [B, (A, C) - (C, A)] „ (210) 

which is a consequence of the Jacobi identity on the Lie bracket [•, •] combined with 
properties (126) and (127) of the bracket (•, •) := [ad^)-, ■], holding for m-valued 
functions A, B, C. To proceed, it is convenient to split up the local and nonlocal 
parts of H on (202). The nonlocal part consists of, via (126), 

-ad{e x y 1 [ad{e x )v 1 D x 1 [ad{e x )v,ad{e x y 1 w { f > ]\\] = [v, D~ x [v, (211) 

which yields three terms: First, 

[^>- D x 1 [ t, '^x]||] =aA{v)[v,v x ]\\ (212) 

follows from the elementary identity D x [v,v x ] — [v,v xx ]. Second, 

[v.D^lv, [aA{e x )v,v x ]±\\\\ = ad^D' 1 ^, {v,v x )]\\ = ^ad(u)[u, (u,u)]|| (213) 

is obtained from the Lie bracket relation (67) followed by the identity 
D x [v,(v,v)}\\ =3[v,(v,v x )]\\ arising from (210). Third, 

-i[t,,£;V[ad(e>,[ad(e>,t;]||]]||] = - |ad(w)Z> J 1 [(w, v) „ , (u, w)] „ =0 (214) 

arises from the same Lie bracket relation combined with (210). Next the local parts 
of H on (202) are straightforwardly given by expanding out 

X~\D x wf) (215) 
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and 

-ad(e :c )- 1 [ad(e ;c ) U ,ad(e :c )- 1 n7( 1) ]. (216) 

Then (205) directly follows by combining these terms with the previous ones. 
The derivation of (206) is similar. 

Now the +1 flow is given by v t = + which is equivalent to the evolution 
equation (187) with h± = and zuj_ = zu^\ This equation can be written out 
in a fairly short explicit form by means of the brackets [•, •] and (■, ■) together with 
the following linear maps 

ad(u)_L := [v,-]±, ad" L (w) := ad^)" 1 [ad( e;I >, (217) 

defined on the vector space mi ~ f)j_. 
Proposition 5.6: 

v t = v x + X~ 1 {v xxx ) — X~ l {^{v, (u, w)|| ) — [v x , v)±) x — (X~ 1 (v xx ), v)± 

-ad(v)±[v,v x + ^(v,v)±]\\ + ^ad ± (v)[v, {v,v)\\] - a,d ± (v) 2 v x (218) 

is the +1 flow in the hierarchy (193). Up to a convective term, this evolution 
equation describes a ii| -invariant mKdV soliton equation on the m±-valued flow 
variable v, with two compatible Hamiltonian structures 

v t -v x = H{5?)/5v) = £(5<E/5v) (219) 

given by 

& = J {-\\v x ? + \{v x , (v,v)) + ±\(v,v)\\\ 2 )dx, <S = - J{v,X{v))dx (220) 
where 

£ :='R'H = -{&d{e x )- 1 K,) 2 . (221) 

The mKdV flow equation (218) corresponds to a geometrical motion of the curve 
"f(t,x) obtained through the identifications (71), (183), (184) arising from m <-> 
T 7 M as defined via a ii-parallel moving frame along j(t, x). Thus, in terms of the 
principal normal variable (183), the mKdV flow is given by 

h± = X- 1 {v x ) = X-\V x v-[u,v])=X- 1 V x v-X- 1 [u,X- 1 [u,e x }] 

= X~ lr D x v - (X^ 1 &d(v)) 2 e x (222) 

and 

H = \[X-\v),X-\ & &{e x )v))\ ] =-\[X-\v),[X- 1 {v),e x ]\ 

= -±( & &(X-\v)fe x \ (223) 

where, recall, v <-> N = V X 7 X , V x v <-> V X N = V 2 j x , e x <-> X — j x . Also note 

X -&d x (X) 2 = -ad( 7;E ) 2 := X 1 (224) 
which represents an invertiblc linear map on T 1 M. Hence 

h ± + h l{ =e t ^ X^V^ - (X^ad^fX -^(^(X^NfX^ = 7t (225) 
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yields the following G-invariant geometrical equation of motion 

7t = X^Vlix - (A^ad^V^))^ - i(ad x (Ar 7 _1 V x 7x) 2 7x)||, It^ Is- = 1 

(226) 

which will be called a non- stretching mKdV map on M = G/H . 

Remark: It is important to emphasize that the mKdV evolution equation (218) 
and the associated geometric map equation (226) are bi-Hamiltonian integrable 
systems that describe multicomponent generalizations of the scalar mKdV equation. 
In particular, the evolution equation (218) possesses a hierarchy of higher-order 
commuting symmetries and higher-order conserved densities corresponding to the 
vector fields (194) and Hamiltonians (199) for n > 1. 

Both the mKdV map (226) and evolution equation (218) greatly simplify in the 
case when the linear map X — — ad^) 2 is a multiple of the identity, 

X = xidj_, X = const (227) 

on the vector spaces mj_ — f)x- From Proposition 2.1, note that (227) implies the 
Lie bracket relations 

[f)x,f)x] C J)||, [m_L,mj_] C f)||, [f)i,mi]Cm||, (228) 

so consequently 

[;-]± = (;-)± =0, ad(«)±- = ad^)- = 0, onmi~f)i, (229) 

while 

ad^v)- = X- 1 Sid(X(v)) ± - = ad(u)- , on my. (230) 
This implies, in addition, the identity 

(A, (B, (J)) - (C, (B, A)) = X [B, [C, A]] (231) 

holding for functions A, B, C with values in mi ~ f)j_. 

As a result, when property (227) holds, the mKdV evolution equation (218) can 
be shown to reduce to 

vt - XVx = v xxx - ^{v x , (v, v)) = H(SH W /5v) = £(SH (Q) /Sv) (232) 

with a factor x absorbed into the time derivative, and = — ^\v\ 2 , = 

— ^\ v x\ 2 + ||( W 7 W )| 2 ! where the bi-Hamiltonian structure becomes 

n = X - 1 D x -[v,D- 1 [v,-}}, j = D x -(v,D- 1 (v,-)), £ = HJH, (233) 
and where the recursion operator is given by 

n = nj = x" 1 (d 2 x - x[v, [v, ■]]-(«,(«,- )) + x[v,d- 1 [v x ,-}} 

-{v x ,D~\v, ■ )) + [v, D-\v, (v, D-\v, ■ ))]]) (234) 

(while the adjoint operator 1Z* is given by simply switching the brackets -^/xl-,-] 
and (•, •) in TV). 

The mKdV map equation (226) similarly simplifies to 

It = V 2 7x - ^X _1 ad x (V x 7x) 2 7x. (235) 

These results (232)-(235) encompass the two versions of mKdV vector evolution 
equations and curve flows in M = SO(N + l)/SO(N),SU(N)/SO{N) derived 
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in [3], as well as the their complex- valued generalizations obtained from M — 
SO(N + l),SU(N) in [4]. In particular the scalar mKdV equation and its bi- 
Hamiltonian structure is given by (232)-(235) in the case M = S 2 . 

5.3. Remarks on higher-order flows. Compared to the mKdV flow (218), the 
+n flow in the hierarchy (193) for n > 2 is a mKdV evolution equation of higher 
order 2n + 1 on the flow variable v. The corresponding curve flows on j(t, x) 
in M = G/H are geometrically described by higher-order versions of the non- 
stretching mKdV map equation (226). All of these flows and associated geometric 
map equations are bi-Hamiltonian integrable systems given by a group-invariant 
multicomponent generalization of the scalar mKdV soliton hierarchy coming from 
curve flows in two-dimensional constant curvature spaces, M = S 2 (cf section 1). 
Like in the scalar case, each flow (193) turns out to be an explicit local polynomial 
expression in v and its x derivatives, for which a general proof should be possible 
by applying the locality results established in [32]. 

There is a direct differential geometric interpretation of the flows (193), stemming 
from the geometrical meaning of the flow variable v = X~ l {e\N) given by the 
principal normal vector N — V ' x "f x along the curve 7 in an _ff-parallcl moving 
frame. Because of the if -parallel property (cf Proposition 3.6), the components of 
v represent differential covariants of 7 relative to the equivalence group HZ C H* 
of the framing. More precisely, these components arc invariantly determined by 
the curve 7 up the covariant action of the (x-independent) group H*^ in the vector 
space C m ~g T x M containing v. In contrast, for the case M — S 2 , there is 
essentially a unique framing and the equivalence group HZ comprises only discrete 
reflections about the unit tangent vector 7^ in the tangent plane T X M ~ M 2 . In 
this case v represents a differential invariant of 7, namely the classical curvature 
invariant. The results in Theorem 5.5 thus show that, in a generalization from 
M = S 2 to symmetric spaces M = G/H, differential covariants of the curve 7 
give a proper setting for a flow variable that geometrically encodes a natural bi- 
Hamiltonian structure. 

5.4. SG flows. In addition to all the local polynomial flows (193) in the hierarchy 
given in Theorem 5.5, there is a nonlocal, non-polynomial flow arising from the 
kernel of the recursion operator (171) as follows. Consider the vector field h±-d/dv 
defined by 

= w± = J(h±) (236) 

which thus gets mapped into TZ(h_\_) = by the recursion operator. This will be 
called the —1 flow, with h± assigned zero scaling weight under the mKdV scaling 
group x — > Xx, v — > \~ 1 v. Geometrically, from relation (185), w± = ejV t 7 2 = 
immediately determines the underlying curve flow on 7 to be the G-invariant 
equation of motion 

V t 7* = 0, |7x| fl = l (237) 
which is recognized as being a non- stretching wave map onM = G/H , with t and 
x interpreted as light cone coordinates for the curve flow. 

The corresponding evolution equation (187) induced on the flow variable v is 
given by h± = hj_ ^ such that zu = J(h± ^) = 0, so hence v satisfies the —1 flow 
equation 

v t = h { ^\ D x h ( - 1] + [ad(e x )v,h ( - 1) -D- 1 [ a d(e x )v,h { - 1 \} ± = (238) 



38 



STEPHEN C. ANCO 



describing a nonlocal flow. The structure of this equation looks simpler if we work 
in terms of the flow-vector variable h± + h\\ — e* = e\jt, and the x-connection 
variable u — ad(e x )v = uj x = oj\-f X} which obey the frame structure equations 
(68)-(70) in Lemma 4.2. Write 



:= h { - 1} + h\- 1} (239) 



and note 



VU ± +ZU 11 =LU t =Lu\~/t = (240) 

where zu 1 - — — a,d(e x )zv± — and nj" = — D x 1 [u, zo 1 -]^ = due to zo± = 0. 

Proposition 5.7: The —1 flow equation (238) is equivalent to the nonlocal 
evolution equation 

u t = &A(e x )h^ (241) 

with 

D x h(~V =-[u,h(~V] (242) 
whose geometrical content expresses the vanishing of the connection (24-0) in the 
flow direction jt ■ 

This flow (241)-(242) possesses the conservation law 

= D x \h^\\ := -(h^\h^) - (/^UP) = \ lt \l, (243) 

geometrically corresponding to 

= V x | 7 tl2 (244) 
admitted by the wave map equation (237). Consequently, a conformal scaling of t 
can be used to put 

l = |fc(-U| = | 7t | fl (245) 

whence the —1 flow is equivalent to a flow with unit speed. 

In the examples M = SO(N + l)/SO(N),SU(N)/SO(N) and M = SO(N + 
1),SU(N), the conformally scaled —1 flow given by equations (241)-(245) was 
shown [3, 4] to admit an algebraic reduction yielding a hyperbolic equation of local 
non-polynomial form v tx = f(v,v t ) on the flow variable v — ad(e;!;) _1 u. Here a 
related hyperbolic equation will be derived that directly generalizes the SG soliton 
equation from the scalar case M = S 2 to M = G/H, without any reductions being 
introduced. 

To proceed, observe that equation (242) has the geometric formulation 

V x h^ = (246) 

where 

V X = D X + [u, ■ } (247) 
is a covariant derivative operator (cf Proposition 4.4) associated with the x- 
connection in the flow. Now consider a local gauge transformation on the x- 
connection, given by a smooth function q : 7 C M — G/H — > H* ~ Ad(H), 
such that the operator (247) is flattened, V x — > qV x q~ x = D x . This has the 
geometrical meaning that the transformed cc-connection vanishes 

u -> quq' 1 - q x q^ := u = (248) 

giving a linear differential equation q x = qu that determines q in terms of u. For 
any such q, it then follows that the transformed flow-vector is constant 

h { ~ 1} -» qh^q- 1 := h, D x h = (249) 



SOLITON EQUATIONS AND BI-HAMILTONIAN CURVE FLOWS 



39 



where h is some constant m- valued vector with unit norm due to property (245). 
As a consequence, equation (246) is solved by 

u = q- 1 q Xl h { ~ V> =q- 1 hq (250) 

in terms of a IP-valued potential q. Substitution of these expressions into the —1 
flow equation (241) now leads to a main result. 

Theorem 5.8: The —1 flow on v = ad(e x ) _1 u = &d(e x )~ 1 (q~ 1 q x ) is given by 
a hyperbolic equation 

qtx ~ q t q~ Y q x = q[e x , q^hq] = [ge x g _1 , % (251) 

for the H* -valued potential q, with D x h — and \h\ m = 1. 

This flow depends on the pair of constant unit vectors e x ,h in m. Without loss 
of generality under a rigid (x- independent) transformation h — > qhq^ 1 if necessary 
(which preserves (250)), the vector h can be taken to belong to the same Weyl 
chamber a*(m) that contains the vector e x in a fixed maximal abclian subspace 
a C m, whereby the flow (251) is parametrized by unit-norm elements in the Weyl 
chamber a*(m) in m (cf Proposition 3.2). Note the invariance group of the flow will 
be iff h = e x . 

A natural Hamiltonian structure can be derived for the hyperbolic flow equation 
(251) starting from the scaling formula (199) extended to n = —1. Let 

H { - 1) := (e x , h^) = {e x ,h { ~ 1] ) = (e^q^hq)- (252) 

Its variational derivative with respect to q is defined by SH^ 1 ^ = 
tr{8q{8H { ~ 1] /8q)) = -(q^Sq, (£i2" (_1) / '5q)q) , which yields 

-{5H ( - 1) /5q)q= [e^q^hq] = &d{e x )h ( - 1] . (253) 

Theorem 5.9: The hyperbolic flow equation (251) has the Hamiltonian structure 

q-\qtx ~ qtq-'q^q- 1 = SH^/Sq (254) 

or equivalently 

qtq- 1 = -D- 1 (q(5H { - 1) /5q)) (255) 

as an evolutionary flow, where D^ 1 is a Hamiltonian operator on the x-jet space 
of ^-valued functions. 

This Hamiltonian structure has an equivalent formulation in terms of the flow 
variable v — &d(e x )~ 1 u by means of the variational relations 

6i2" (_1) = -(Su,5H { - 1) /5u) = -(Sv^H^/Sv), (256) 
Su = ad{e x )Sv = Siq-iq*) = V^q^Sq), (257) 

combined with the variational derivative (253), where 5H 1 " ^ /Su and SH^ ^/Sv 
are t)j_- and m_i_-valued respectively. Note, since (<5w)|| = 0, the relation (257) 
implies (q~ 1 Sq)\\ = — D" 1 ^, (q~ 1 8q)±] \\ and hence 

Su = D^q^SqU + [u, {q^Sq^ - D^[u, (g" 1 ^) = /C^" 1 ^ (258) 

where JC is the operator (74). 

Proposition 5.10: The —1 flow equation (24-1) -(24-2) has the Hamiltonian 
structure 

u t =n(SH i ~ 1) /Su) (259) 
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where ^ = (e x ,q 1 hq) is a nonlocal function of u determined by q x = qu. 
Equivalently, the — 1 flow equation (238) on v — ad(e x ) _1 w is given by 

v t =H(5H { - 1) /Sv) (260) 

with q expressed in terms of v through q x = q[e x ,v]. 

As a consequence, the —1 flow on v thus defines a nonlocal Hamiltonian vector 
field 

h { ~ 1] := {q~ 1 hq) ± = H{5H^~ V> /Sv) (261) 

and likewise 

hf_ 1} := ad(e x )/ li 7 1) = [e x , q^hq] = HiSH^/Su) (262) 

is a corresponding nonlocal Hamiltonian vector field given by the —1 flow on u. 

Remarks: The group-valued hyperbolic equation (254) provides a multicom- 
ponent generalization of the scalar SG equation along with its Hamiltonian struc- 
ture, arising from curve flows (237) in M = G/H. In particular, for M = S 2 
if H* = 50(2) ~ U(l) is viewed as the group of abelian unitary matrices then 
q = exp(i0) can be identified as a {/(l)-valued potential satisfying the SG equation 
6 tx — sin(# — 6) with the Hamiltonian ^ = i3(exp(i#) exp(— id)), 9 =constant. 
From Theorems 5.8 and 5.9 and Proposition 5.10, the generalization of this soliton 
equation in the multicomponent (nonabelian) case describes an integrable SG sys- 
tem for which the associated —1 flow in nonlocal evolutionary form (255) possesses 
a hierarchy of higher-order commuting symmetries given by the vector fields (194) 
and higher-order conserved densities given by the Hamiltonians (199) for n > 0. 
Moreover, these hierarchies can be expressed directly in terms of the potential q 
via the relation v = &d(e x )~ 1 u — ad(e x ) _1 (q~ 1 q x ). 

6. Concluding Remarks 

The main results in this paper on bi-Hamiltonian hierarchies of geometric curve 
flows can be directly extended to flat Riemannian symmetric spaces m = g/t) ~ R™ 
viewed as infinitesimal Klein geometries. One significant geometrical difference will 
be that both the curvature and torsion of the Riemannian connection V will vanish 
in the Cartan structure equations for the framing of the curve flows ~f{t,x) in R™. 
This change does not affect the encoding of bi-Hamiltonian operators, which will be 
the same Hamiltonian cosymplectic and symplectic operators as derived in curved 
symmetric spaces M = G/H, but it does lead to a simpler recursion structure 
for the resulting integrable hierarchy of group-invariant flow equations induced on 
the moving frame components of the principal normal vector along the curves. In 
particular, the higher-order mKdV flows in each hierarchy for M = G/H (when 
V has non-zero curvature) are linear combinations of the ones arising in m = g/f) 
(when V is flat). 

A further extension of the main results in the present paper arises for Riemannian 
symmetric spaces M = G/H ~ K given by any compact semisimple Lie group K, 
where G = K x K is a product group and H ~ K is a diagonal subgroup. Since M 
here is a group manifold, there turns out to be an algebraic Hamiltonian operator 
that is compatible with the cosymplectic and symplectic operators encoded in the 
Cartan structure equations for the framing of curve flows 7(i, x). This leads to 
an enlarged bi-Hamiltonian hierarchy of geometric curve flows, containing group- 
invariant NLS flows in addition to mKdV flows. The same result extends to the 
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flat spaces m = g/l) ~ t for any compact semisimple Lie algebra t (viewed as an 
infinitesimal Klein geometry). Details of this will presented in a forthcoming paper. 

Enlarged bi-Hamiltonian hierarchies of geometric curve flows will likewise arise 
for other symmetric spaces that carry extra algebraic structure, such as hermit- 
ian spaces or quaternion spaces. All such spaces are well-known from Cartan's 
classification [19]. 

Earlier work [3, 4] studying the four symmetric spaces SO(N + l)/SO(N), 
SU(N)/SO(N), SO(N + 1), SU(N) geometrically accounted for all examples of 
0(N — l)-invariant and U(N — l)-invariant multicomponent mKdV and SG equa- 
tions that have been found in symmetry-intcgrability classifications of vector evolu- 
tion equations [34] and vector hyperbolic equations [6]. Building on these examples, 
a subsequent paper will be devoted to a full classification of further examples of 
group-invariant soliton equations derived from geometric curve flows in all classical 
(non-exceptional) symmetric spaces. 

Finally, in light of work done in [10, 11] deriving numerous scalar soliton equa- 
tions from geometric curve flows in planar Klein geometries G/H ~ M 2 , a natural 
interesting question is whether the present analysis using moving parallel frames 
for curve flows in Riemannian symmetric spaces can be broadly extended to gen- 
eral homogeneous spaces M = G/H. Such a generalization would encompass non- 
Riemannian Klein geometries and, accordingly, should yield group-invariant soliton 
equations (along with their bi-Hamiltonian structure) of other than mKdV and SG 
type. 
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